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ABSTRACT 



This thesis report discusses the design, construction and 
calibration of two force-torque transducers for use in a force 
control override of a rate control system. Pre-loaded force 
sensing resistors were used in a computer model to determine 
the number and location of sensors necessary to resolve three 
forces and three moments. An analysis was conducted on this 
full order model to determine redundancy limits. A reduced 
order model was then used to determine the sensor 
configuration required to resolve three forces and only one 
moment. Prototypes of the reduced order model were then built 
in two different sizes, and used to sense and display applied 
forces and moments. 
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I. 



INTRODUCTION 



The objective of this thesis is to design, construct and 
calibrate a pair of reduced order force-torque transducers 
based on force sensing resistor technology. This is one 
branch of a two branch robotics research project. The work 
done in parallel is a separate study of how a hybrid force 
override rate control system can be used to enable operation 
of robotic manipulators. 

The motivation for this research is that robotic 
applications such as the Space Shuttle Manipulator Arm require 
more than just rate and position control to perform certain 
assembly tasks. The ability to detect and control interaction 
forces and moments is necessary in poorly defined environments 
such as space. To accomplish various tasks without damage to 
the surrounding objects or the manipulator, the contact forces 
must be sensed and used in a feedback loop to regulate the 
conventional rate and position control. Further, with the use 
of force override control systems, human intervention, which 
may be hazardous, impractical or even impossible, could be 
entirely eliminated. This offers the potential for reducing 
costs, protecting human life and enabling the accomplishment 
of tasks never before possible. 
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In order to develop a practical, simple and durable 
system, the design of the force-torque transducer in this 
project is based on the Force Sensing Resistor™ (FSR) rather 
than the strain gauge. FSRs provide a relatively inexpensive, 
vibration insensitive, durable and adaptable source of 
resistance changes from various force applications. To harness 
this technology in the design of a force-torque transducer, it 
was necessary to determine a basic configuration which enables 
the desired forces and moments to be extracted from various 
sensor signals. A design which enables the resolution of three 
forces and three moments was developed. This design was then 
modified to provide three forces and only one moment in the 
likely scenario of applications requiring only these four 
commonly used degrees of freedom. Also, the redundancy of the 
basic six degrees of freedom transducer was analyzed keeping 
in mind the practical use of such a device. A scheme allowing 
multiple levels of redundancy, in the event of individual 
sensor failure, was developed to demonstrate the functional 
reliability of a force-torque transducer designed using force 
sensing resistor technology. 

Finally, the four degree of freedom transducer and 
associated electronic interfaces were constructed and tested 
in two versions. One was a design requiring no alterations to 
the readily available FSR. The other required a new, shorter 
FSR not currently available from the manufacturer. To overcome 
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this obstacle the available FSRs were modified in order to 



build a force-torque transducer of reduced proportions. This 
was done due to the payload limitations of the PUMA 560 
robotic manipulator arm. 
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II 



PRELIMINARY WORK 



A. FORCE OVERRIDE RATE CONTROL 

1 . Concept 

The general concept of a force override rate control 
system is that the contact forces between a robotic 
manipulator end effector and the object being moved can be 
used to generate a feedback signal. This signal is in the 
form of a voltage change which is used to alter the control 
input signal of a conventional rate or position control 
system. This gives conventionally controlled manipulators the 
ability to reflect the forces applied by the operator on the 
control device. 

2 . Background 

This concept has been tested using a single degree of 
freedom hydraulic system [Ref.l:pp. 5-12], a single degree of 
freedom PUMA 560 manipulator system and a three degree of 
freedom PUMA 560 manipulator system [Ref. 2 ;pp. -27-54 ] . Figure 
1 provides a general description of the hydraulic system [Ref. 
l:p.6]. Figure 2 shows the general arrangement of the PUMA 
system [Ref. 3:p. 23]. These systems utilize resistance strain 
gauges to generate voltage signals which are proportional to 
the end effector contact forces. 
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Figure l 

Hydraulic Force Override Rate Control System [Ref. 2 :p. 28] 




Figure 2 

One DOF Force Override Rate Control of a PUMA 560 
Manipulator [Ref. 2:p. 34] 
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In general a strain gauge is used to determine the 
actual stress experienced on the surface of the member to 
which it is bonded. The device consists of a folded wire 
assembly that exhibits a resistance change as the length and 
cross sectional area change according to: 

R=pL/A (1) 

Where p is the resistivity of the material. L and A represent 
the wire length and cross sectional area. Figure 3 shows an 
enlarged view of a strain gage. These resistance changes are 
small and temperature effects must be compensated for by use 
of a second unstrained gauge as part of the system. Also, the 




gauge must be bonded to the stressed member in such a way that 
all the strain is transferred to and detected by the gauge. 
Strain gauge rosettes consisting of two or more strain gauges 
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with different orientations are used when stresses are 



expected in more than one direction. Hooke's law is used to 
calculate the stresses from the strains which are proportional 
to the applied loads. 

Force-torque transducers based on strain gauge 
technology have been applied in the design of wrist force 
transducers. However, these tend to be expensive and complex 
requiring the part on which the strain gauge is mounted to be 
capable of deflection [Ref. l:p. 15]. Also, such systems have 
a high susceptibility to noise since the small resistance 
changes require very high amplifier gains in order to produce 
a useful voltage signal. Hence, an alternative source of 
voltage output proportional to force input was sought which 
was low cost, robust and simple to construct, calibrate and 
use . 

B. SENSOR DEVELOPMENT 

l. Force Sensing Resistor 
a. Description 

(1) Construction 

A Force Sensing Resistor™ is a device 
manufactured by Interlink Electronics of Santa Barbara, CA. 
Figure 4 shows the construction of a typical FSR. A sheet of 
polymer with a layer of semiconductive ink is faced with 
another sheet of polymer overlaid with a conducting pattern 
consisting of a set of interconnecting electrodes. These 
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Figure 4 

Typical Force Sensing Resistor [Ref. 4:p. 1] 



electrodes are shunted by the semiconducting ink when the 
layers are brought together . [Ref . 4:pp. 1-2] 

(2) Characteristics 

The FSR has a very high no load resistance on 
the order of 1 Mohm or more. As the layered device is 
compressed the resistance drops proportionally. Figure 5 shows 
a typical plot of resistance verses force which is 
significantly more sensitive than the output of a strain 
gauge. Thus, the susceptibility to noise is immensely 
reduced. [Ref . 4:pp. 1-2] 

Figure 6 displays the results of a longevity 
test conducted by the manufacturer demonstrating the 
durability of a typical FSR. Interlink Electronics also 
reports that the standard temperature range is up to about 
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170°C. Sensors capable of withstanding higher temperatures 
can be manufactured. FSRs are also relatively insensitive to 
humidity . [Ref . 4:p. 3] 



FSR Force'Resistance Characteristic 




Figure 5 

Typical FSR Force/Resistance Characteristic [Ref.4:p. 1] 



FSR Longevity Test 




Test Conomoos: 1/2* cUa. circular prooe 
10 million stnkes ai 52 
0.5 sec/reoetmon 



Figure 6 

FSR Longevity Test [Ref.4:p. 3] 
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FSRs can be readily manufactured in various 
sizes as single sensors or arranged in arrays. It is also 
possible to vary the force range and the resistance range to 
meet specific requirements during fabrication. Interlink 
Electronics suggests that it is best to use the mechanical 
design to establish the force within a useable range. [Ref. 
4 : p. 2 ] 

An FSR's performance is somewhere between a 
true force sensor and a pressure sensor. For a given force, 
a pressure sensor yields an output inversely proportional to 
the area of the applied force. Whereas a force sensor will 
yield a constant output regardless of the area and 
distribution over which the given force is applied.. If the 
force distribution area is smaller than the FSR active area 
but large compared to the spacing between the interconnecting 
electrodes, then the resistance output will vary approximately 
as the reciprocal of the square root of the applied force 
area. However, when the area of the applied load is larger 
than the FSR active area the device can - be used as a 
pressure sensor. [Ref. 4:p. 2] 

Sensor electrical interfaces can be simple 
because typical resistance changes are relatively large; 
therefore, these circuits do not require a resistance bridge 
as do strain gauge circuits. Also, the impedance is almost 
entirely resistive. Depending on the circuit used, sensor 
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response signals are typically on the order of 0-10 
volts . [ Ref . 4 : p . 4 ] 

Force sensing resistor technology offers many 
advantages over conventional strain gauge technology. In an 
attempt to develop the force override rate control system as 
a low cost, rugged and credible concept, a suitable force- 
torque transducer must be developed. Hence, the 
implementation of readily available FSRs is a logical choice 
for a simple, robust and inexpensive transducer design. 

2 . Transducer 
a. Principle 

The use of force sensing resistor technology was 
suggested to the project sponsor at NASA. Figure 7 shows the 
sponsor's suggested conceptual design for a suitable force- 
torque transducer using readily available FSRs [Ref. l:p. 16]. 
Initially, it was proposed that the sensors be mounted on an 
inner cube which would be surrounded with a semi-compliant 
material like RTV and encased in a larger cube. When the 
joystick is moved, some of the FSRs will experience local 
normal forces and thus a change in resistance. When the 
appropriate electronic interface is used, a corresponding 
voltage change will be produced resulting in sensor signals 
which reflect the applied forces and moments. Ultimately, the 
sensors would be located in such a way that the applied forces 
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and moments could be detected and decoupled for use in a force 
override rate control system. (Ref. 3:p. 11] 




Figure 7 

FSR- Based Transducer Conceptual Design [Ref. l:p. 17] 
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b. Sensor Equations 



The original mathematically generalized description 
of the problem was: 




a l,2 a l,3 * 

• • • 

• • • 

• • • 

• • • 

a n,2 a n,3 * 



a l,12 













( 2 ) 



Vi 2 



My + 

My- 

M z+ 

M, 



In this equation the sensor vector, S, represents the voltage 
output for each sensor when the force and moment vector, F, is 
applied, and n represents the number of sensors used which 
depends on the selected configuration. The positive and 
negative directions for each force and moment component had to 
be accounted for because the sensors were only used to detect 
applied compressive loads. Hence, the objective was to 
determine the right number and location of sensors and the 
associated A matrix in order to resolve the 12x1 force and 
moment vector using: 

F=A' 1 S (3) 
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The elements of the A matrix would have to be determined using 
a calibration procedure which measured each sensor signal 
after applying known forces and moments individually. [Ref. 
l:p. 18] 

Since the number of sensors could be more than 
twelve, the A matrix may not be invertible. In which case the 
pseudoinverse must be utilized according to: 

F= (A t A) 'A t S (4) 

If the square matrix A T A is of full rank (i.e. a rank equal to 
the length of the force and moment vector) , it may be inverted 
and used as shown in Equation 4 to identify the forces and 
moments. [Ref. l:p. 19] 

The initial model of a control device to detect the 
forces and moments was that of a planar joystick shown in 
Figure 8. An equation was developed for each sensor using the 
principle of superposition and modeling each FSR as a linear 
spring. The principle of superposition was used because the 
application of a force, P, at the top of the joystick would 
cause the device to translate and rotate. Hence, the 
following equations account for the contribution of both 
effects to the normal force on each spring: 



F,=P/4(l+d/a) 


(5) 


F 2 =P/4(l-d/a) 


(6) 


F 3 =-P/4 (1+d/a) 


(7) 


F 4 =-P/4 ( 1-d/a) 


(8) 
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Figure 8 

Planar Joystick [Ref. 3:p.20] 



By applying a pure moment about the y axis the equations 
become: 

F,=F 4 =M/4a (9) 

F 2 =F 3 =-M/4a (10) 

The sensor coefficient matrix A can be determined for a 
variety of three dimensional configurations by extending this 
methodology to a cube and determining the effects of each 
applied force and moment. (Ref. l:pp. 17-21] 
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c. Sensor Placement 



Since the A matrix is dependent upon the location 
of the sensors on the cube, the first task in this iterative 
design process was to pick a number and distribution of FSRs 
on the cube. Then by applying the forces and moments of 
interest, the components of the A matrix were determined using 
the principle of superposition outlined earlier. The A T A 
matrix was then computed and its rank checked. If the results 
indicate a rank equal to the length of the force and moment 
vector, then that sensor configuration was considered capable 
of being used to resolve the forces and moments. If the rank 
was greater than the vector length, then too many sensors were 
being considered. If the rank was less than the vector 
length, then too few sensors were being considered. [Ref. 3:p. 
14] 

d. Early Results 

The iterative process for determining the right 
number and location of sensors was carried out for a full 
order transducer consisting of 6 forces and 6 ^moments. This 
accounted for three dimensional forces and moments with 
sensors only capable of registering a compressive load. A 
possible placement pattern which resulted in an A r A matrix 
with a rank of 12 is shown in Figure 9. [Ref. l:pp. 29-30] 

The iterative process was then carried out in a 
similar fashion for a reduced order transducer capable of 
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resolving the three dimensional forces as before, but a moment 
only about the z axis. This was done because not all robotic 




Figure 9 

Arrangement for Full Order Transducer [Ref. l:p. 30] 



manipulator tasks require moments in all three directions. It 
is common to only require a torque about a single axis in 
addition to the forces along all three directions. Again, 
since the FSRs are only responsive to compressive loads, a 
separate sensor had to be allowed for the negative and 
positive directions of each of these components. The 
potential sensor layout with a rank of 8 is shown in Figure 
10. [Ref. 1 :pp. 31-33] 
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It is apparent that in order to reduce the number 
of sensors required, a design modification is necessary which 
addresses the response of FSRs only to compressive loads. 




Figure 10 

Arrangement for Reduced Order Transducer [Ref. l:p.33] 



This could, in effect, halve the force and moment vector 
length by not requiring sensors for negative direction forces 
and moments . 
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III. PRE-LOADED SENSOR TRANSDUCER DEVELOPMENT 

A. SINGLE FSR PROTOTYPE 
1* Objective 

In an effort to resolve the problem of the FSR's 
inability to register a signal when a tensile load is applied, 
the concept of pre-loading the sensor was employed. This 
required the sensor to be loaded compressively with some 
initial force. The change in resistance caused by adding or 
removing normal force applied to the FSR is used to generate 
a proportional change in voltage. Provided that the initial 
pre-load is sufficient, a useable range is created in which a 
response is generated to both tensile and compressive force 
application. 

2. Design and Construction 

In order to test this concept a device had to be 
created which would show a reasonable range of resistance 
change for varying compressive and tensile- normal force 
application. 

a. Mechanical Device 

Figure 11 shows the prototype used to test the pre- 
load concept. The device was constructed using standard .25" 
aluminum plate. As shown .25" of neoprene material was used 
to transmit the applied load to the FSR which was mounted 
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using Scotch™ 467MP High Performance adhesive. The FSR was 
.25" square, product #30-301, readily available from Interlink 
Electronics™. A standard 1.5" 10/32 hex-head machine bolt 
with a locking nut was to apply the initial pre-load. A 
detailed drawing is provided in Appendix A. 




Figure ll 

Single Pre-loaded FSR Prototype 
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b. Electronic Interface 



In order to obtain a voltage change proportional to 
the applied load, the electronic interface pictured in Figure 
12 was used. The circuit consisted of a +10 volt reference 
source, various values of and an LM324 op-amp. The output 
voltage can be written as: 

^out = ^r£f/ 2*[1+R 0 /R fsr ] (11) 
The expected voltage swing is from V^/2 to V,^. [Ref. 5:p. 8] 



R G 




Figure 12 

FSR Current-to-Voltage Converter [Ref. 5:p. 8j 
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3. Response Test 

a. Data Collection 

The prototype was configured to sit on a table top 
while a compressive load was applied by adding weights to a 
tray on top of the pre-load set bolt. Tensile loads were 
applied by turning the device over and suspending it by 
supporting the base on two cross members across a sufficient 
gap between two adjacent and parallel tables. The weights 
were then added to a tray suspended from the pre-load set 
bolt. 

Initially, the masses, ranging from 4.17 to -1.97 
lbm, were measured in English units and the corresponding 
voltages were read using the OMEGA BENCH™ data acquisition 
system and APPLE MACINTOSH™ microcomputer configured as a 
voltmeter. The first test used 4.7 kohms for the value of Re 
in an effort to produce a curve similar to that shown in 
Figure 13. This was the response for a standard FSR using a 
similar interface obtained by Interlink Electronics™. [Ref. 
5:p. 7] 

In order to more fully understand the nature of 
pre-loaded FSRs, similar tests were conducted using the 
following values for R*, in kohms: 2.5,5.0,7.5,10.0,12.5. The 
masses used during these tests were measured in metric units 
for convenience. 
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Rc Values 



- — 75k 

< 5 — 4.7k 

* — 15k 



F vs. V for part # 174 
Interlink Force Tester 
1 cm 2 circular flat 
metal actuator. 

0 200 400 600 800 1000 

Force (g) 

Figure 13 

Force vs. V oirr Curves [Ref. 5:p. 7] 



b. Results 

The data taken and the corresponding plots of 
applied force versus voltage response are provided in Appendix 
A. The plot resulting from the first test is shown in Figure 
14 . The curve resembles the expected results shown in Figure 
13. The plot clearly shows two nearly linear regions 
separated by the initial no load set points. In subsequent 
tests the initial set point was physically unchanged. 
However, different values for Rq resulted in curve shifts. 
This is seen in the variation of the no load voltage levels in 
accordance with the various values for Rq. Also, depending on 
the value of the initial voltage, the useable range of the 
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pre-loaded sensor was affected. Too high a value for Rq 
caused positive saturation as in the 12.5 kohm case. The 
conclusion drawn from these tests was that the FSRs configured 
in a simple pre-loaded manner with a simple electronic 
interface are capable of registering a potentially useable 
voltage change which over limited regions varies almost 
linearly with applied loads in both the compressive and 
tensile directions. 

B. INITIAL TRANSDUCER CONFIGURATION 
1 . Arrangement 

Since each FSR can now be used to detect normal force 
loading in both the compressive direction and the tensile 
direction, the 12x1 force and moment vector in Equation 2 is 
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reduced to a 6x1 vector. The mathematically generalized 
description of the problem is now: 
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Where S represents the sensor voltage output and F is the 
force and moment column vector. The quest is still to 
determine a number and configuration of sensors such that the 
6x1 force and moment vector can be resolved by inverting the 
A matrix according to: 

F=A'*S (13) 



If the sensor configuration is such that only 6 sensors are 
used, then the A matrix will be square and the pseudoinverse 
will not be required as shown previously in Equation 4. 

Figure 15 shows a three dimensional transducer built 
upon the planar joystick concept discussed in the previous 
chapter. The plan was to develop the equations necessary to 
determine the A matrix with all the sensors shown. Then, 
systematically cut out the signals from various groups of 
sensors and check the rank of the corresponding A matrix. If 
the rank was greater than or equal to six, then the 
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configuration was considered capable of resolving the 6 
components of the force and moment vector. The goal was to 
determine a configuration which uses the minimum number of 
sensors and has a corresponding invertible A matrix. 




Figure 15 

Initial Transducer Configuration 
2 . Sensor Equations 

Using the principle of superposition mentioned earlier 
for the planar joystick, equations were developed for forces 
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applied in the x, y, and z directions as well as the moments 
about each axis. These forces and moments were considered to 
have been applied at the top of the joystick handle. First, 
a force applied in the positive x-direction, F x/ results in a 
translation of the inner cube as well as an incidental moment 
about the y axis. Since there are eight sensors on the x 
faces, the translation force contribution is simply: 

F ' =F X / 8 (14) 

The sign of this force equation for each sensor depends on 
whether the sensor is on the positive or the negative x face. 
The incidental moment created by F x causes activation of 16 
sensors, eight on the x faces and eight on the z faces. Hence, 
the resulting applied load due to this incidental moment can 
be determined for each affected sensor according to: 

16F”a=M (15) 



where , 



hence. 



M=F x d 



(16) 



F"=F x d/16a (17) 

Again, the sign of this force equation for each sensor depends 
on its location. Since this force is a result of the 
incidental moment, the location of the sensor with respect to 
the axis of rotation, the y axis in this case, determines the 
sign of F". 
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Finally, with the signs determined correctly for each 
activated sensor's force equations, the normal force 
experienced by each FSR resulting from the application of F x 
is the sum of the translational force and the force from the 
incidental rotation: 

F=F'+F" (18) 

This rationale can be repeated for a force applied in 
the y direction, F y . However, a force applied in the z 
direction, F z , would not result in an a rotation of the inner 
cube. Therefore, the resulting normal forces on the affected 
sensors would only be from the translation of the inner cube 
with no incidental moment. 

For a pure torque applied at the end of the joystick 
in the x direction, M x , there would be sixteen sensors 
affected, 8 on the y face and 8 on the z face. The resulting 
force magnitude for each activated sensor is: 

F=M x /16a (19) 

The sign for each of these activated sensor equations depends 
on its location with respect to the axis of rotation. Pure 
moments applied about the y axis and about the z axis. My and 
M z respectively, result in similar force equations for each 
affected sensor. 

The equations for each sensor activated as a result of 
the application of the individual components of the force and 
moment column vector make up the A matrix in Equation 12. The 
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complete equation for this configuration, showing the 



individual elements of the A matrix, is as follows: 
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3 . Sensor Placement 

This A matrix, resulting from the configuration shown 
in Figure 15, was entered into the software program MATLAB™ 



with the assumed values for d and a as 4.5 and 1.5, 
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respectively with units of length. The the rank of the full 
matrix was then determined to be six. This full rank meant 
that if all sensors were used, the pseudoinverse could be used 
to resolve the components of the force and moment vector, F. 

Next, the sensors on the negative x face were 
eliminated by setting the A matrix elements in corresponding 
rows equal to zero. In this case, it was simulated that 
sensors labelled -1, -2, -3, and -4 had no output. Thus, the 
elements in rows 13 through 16 of matrix A were set equal to 
zero. The resulting matrix had a rank of six. Hence, this 
configuration could be used to resolve the components of F. 

This process was repeated in a similar fashion by 
first eliminating additional sensors from the y face, and then 
from the z face. The result was a configuration with sensors 
only on the positive x, y, and z faces yielding an A matrix 
with a rank of six. Since this would involve a total of 12 
sensors, the A matrix would not be square. However, the 
pseudoinverse could be used to resolve the forces and moments. 

In an effort to reduce the number of sensors required 
even further, this iterative process of removing sensors from 
the configuration was continued. The next removal was that of 
the positive z face sensors. This resulted in an A matrix 
rank of five which meant this sensor configuration could not 
be used. Hence, two sensors were added back to the positive 
z face. These were the diagonally arranged sensors numbered 
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10 and 12 in Figure 15. The resulting A matrix yielded a rank 
of six. This indicated a useable combination. 



diagonal patterns from both of the positive x and y faces. 
The sensors numbered 2, 3, 5 and 8 were eliminated. The 
resulting rank was five indicating an unusable combination. 
Hence, the iterative process was continued to find a useable 
combination of only six sensors. 

4 . Results 

Experimenting with diagonal combinations lead to the 
useable configuration shown in Figure 16 after only three more 
iterations. This is not the only possible combination, but it 



Using this rationale, two sensors were also removed in 
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Figure 16 

Initial Minimum Sensor Transducer Configuration 
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is one that yields a square A matrix which is directly 
invertible to determine the force and moment vector, F. For 
fewer than six sensors, the rules of matrix multiplication 
prohibit the formation of the pseudoinverse. Hence, no fewer 
than six sensors were considered in any configuration. 

C. COMPARISON WITH PRELIMINARY WORK 

1 . Arrangement 

For comparison purposes, an unsuccessful configuration 
analyzed in the preliminary work for a full order transducer 
without using the pre-loaded sensor concept was reexamined 
using pre-loaded sensors. This arrangement is shown in Figure 
17 [Ref. l:p. 25]. The goal was to determine what 
modifications, if any, had to be made to create a 
configuration using only six sensors which would produce an A 
matrix with a full rank. 

2 . Sensor Equations 

Using this configuration and the fact that the 
individual FSRs are capable of measuring tensile and 
compressive loads, a new set of governing force equations were 
developed. These equations were determined using the principle 
of superposition analogous to the methods used in the previous 
section by individual application of each component of the 
force and moment vector, F. 

Since only 12 sensors were considered in this initial 
configuration, in contrast to the 24 of the previous section, 
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Figure 17 

12 Sensor Preliminary Work Configuration [Ref. l:p. 25] 



the complete equation is as follows: 
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Since there are half as many sensors over which the forces are 
distributed, the resulting elements of the A matrix are twice 
as large as in Equation 20. This has absolutely no effect 
since multiplying the A matrix by a scalar will not change its 
rank. 



3. Sensor Placement 



Using MAT LAB™ 


as before, 


the full 


A matrix was 


entered and its rank 


determined 


to 


be six. 


The length 


dimensions, d and a, 


where assumed 


to be 


5.0 and 1.0, 



respectfully. This full rank indicated that the configuration 
could be used to resolve the six components of the force and 
moment vector, P. In an attempt to optimize the number and 
placement of sensors, the sensors on the negative x, y, and z 
faces were eliminated, leaving active sensors numbered 1, 2, 
5, 6, 11 and 12. The resulting configuration yielded an A 

matrix with a rank of only five. Since this is a rank 
deficient condition, another attempt was made by eliminating 
the sensors on only the negative x and y faces, and the 
sensors on the positive z face. This configuration left 
sensors numbered 1, 2, 5, 6, 9 and 10 active. The rank of the 
resulting matrix was five, again a deficient value, indicating 
an unusable combination. 

4. Results 

It was intuitively determined that in order to develop 
a combination utilizing the sensors on only the positive faces 
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or only the negative faces, this configuration may need a 
slight alteration. Hence, it was abandoned in search of a 
similar yet useable configuration. 

D. MODIFIED PRELIMINARY CONFIGURATION 
1. Arrangement 

Considering the previous configuration, Figure 18 
displays a similar arrangement. The only difference is the 




Figure 18 

Modified 12 Sensor Preliminary Work Configuration 

orientation of the z face sensors. Diagonally located sensors 
numbered 11 and 12, located on the positive z face, have been 
rotated 90° clockwise if viewed looking in the positive z 
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direction. Similarly, sensors numbered 9 and 10 have been 
shifted 90° clockwise if viewed in the positive z direction. 

2 . sensor Equations 

In this case the sensor equations which comprise the 
elements of the A matrix were very similar to those of the 
previous section. Since the number of sensors was the same, 
the only change was in the signs of the various force 
contributions as a result of the changed sensor orientation. 
The complete equation describing the new situation is: 
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As expected the only changes were the signs of the elements 
corresponding to the application of F y and M, on sensors 



numbered 9, 10, 11, and 12. 
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3 . Sensor Placement 

MATLAB™ was used to first check the rank of the A 
matrix with all 12 sensors activated. This value was six, 
again indicating a useable configuration. Next, the sensors 
located on the negative x, y, and z faces were again 
eliminated. The rank of the associated matrix was determined 
to be six. With this rank, the configuration is useable. 

4. Results 

A rank of six indicates a full rank allowing the A 
matrix to be directly inverted to determine the forces and 
moments according to Equation 3. Since the minimum number of 
sensors possible is six, this configuration was determined to 
be one of many acceptable choices. Using the trial and error 
methodology it is possible to determine other configurations 
using only six sensors which yield an A matrix capable of 
direct inversion. 
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IV. REDUNDANCY ANALYSIS 



A. GOAL 

Given a transducer with four sensors mounted on each face 
of the inner cube, there is more than one combination of the 
minimum number of sensors, six, which will provide an A matrix 
with a full rank used to resolve the six components of the 
force and moment column vector, F. It would be advantageous 
to develop a system that would automatically find another 
useable combination of remaining sensors necessary to continue 
operation in the event of a sensor failure while in use. 
Hence, the aim was to develop an algorithm which would search 
for and identify various combinations of six sensors that 
could be used to identify the force vector, F. This was to be 
done with a 24 sensor transducer utilizing some initial 
useable configuration of six sensors subject to random FSR 
failures which would drive the system to search for an 
alternate configuration. The number of possible combinations 
would be described by the following equation: 

C(n,r)= nl/(n-r)!r! (23) 

Where n represents the total number of sensors and r 
represents the number of sensors in each group necessary to 
resolve the force and moment vector. [Ref. 6:p. 1-8] 
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In this case, n equals 24 and r equals 6. After doing 
the calculation, the resulting number of combinations is 
134,596. Potentially, this many combinations would have to be 
analyzed to find another configuration which would yield a 
full rank A matrix. However, as sensors fail the number of 
possible combinations decreases. It was decided that with this 
finite limit to the number of possible search iterations, a 
FORTRAN code would be developed to mimic the methodology of 
the previous chapter in finding a useable combination. 

B . THEORY 

1. Sensor Configuration 

As a prelude to developing this FORTRAN code, a 24 
sensor configuration had to be developed with a convenient 
numbering system. Also, it was decided to use a configuration 
that could be initially broken into four separate, six sensor 
configurations which were symmetrical and individually capable 
of resolving the three components of force and three 
components of moment. The purpose of these subsystems was to 
initially provide four layers of redundancy. In the event of 
a sensor failure, a quick switch to another subsystem would 
allow continued use of the joystick without the delay of 
randomly searching for other possible combinations of six 
working sensors from the remaining sensors. This could be 
repeated until no intact subsystems remained. 
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To do this, a 24 sensor transducer, similar to that 
shown in Figure 15, was used. The initial numbering system 
thought to be best was a simple numbering of the sensors from 
#1 through #24. This is shown on the inner cube faces in 
Figure 19. Since it was shown in the previous chapter that 




First Sequential Cube Face Numbering Scheme 

only half the sensors on the positive cube faces were required 
to develop a useable configuration, the symmetrically 
analogous configurations were obvious choices for other 
possible configurations. Hence, the four separate 
configurations were comprised of two positive face systems and 
two negative face systems. As shown in Figure 19, these 
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systems were respectively labeled: A-positive, B-positive, A- 
negative and B-negative. 

To develop a numbering system such that the complete 
24x6 A matrix could be conveniently broken into the four 
separate 6x6 smaller matrices each with a rank of six as a 
result of the four separate six sensor configurations, a trial 
and error procedure was necessary. First, the numbers of the 
sensors were arranged such that A-positive system had sensors 
#l-#6; the A-negative system had sensors #7-#12 : the B- 

positive system had sensors #13-#18; and the B-negative system 
had sensors #19-#24 as shown in Figure 19. Then, in a manner 
similar to that in the previous chapter, the equations 
representing the normal forces on each sensor were developed 
as a result of the individual application to the components of 
F. As in the previous chapter, these equations were placed in 
the A matrix and values for the length dimensions d and a were 
assumed to be 4.5 and 1.5, respectively. The matrix rank was 
determined using MATLAB™ to be six. Next, the four 
subsystems were formed into smaller 6x6 matrices and checked 
for ranks of six. However, the B-positive and the B-negative 
subsystem matrices each had ranks of only five. Hence, this 
configuration was not capable of providing the four levels of 
redundancy initially sought. 

To find a configuration that would yield four useable 
subsystems, the sensor numbering scheme had to be changed. As 



41 



shown in Figure 20, the first such change was to swap #5 with 
#17, #6 with #18, #11 with #23 and #12 with #24. This scheme 
failed to produce the four matrices with ranks of six. In this 




case the A-positive and the A-negative subsystems were rank 
deficient. In an effort to find other arrangements which 
consisted of four useable subsystems sequentially numbered, 
this iterative procedure was continued. Table 1 summarizes 
the results of all the sensor number switches. 

The third scheme, shown in Figure 21, was similar to 
the first scheme except that the switch of #6 with #17 and #12 
with #23 was made. The resulting matrices were of full rank. 
Hence, this was the numbering scheme chosen for the FORTRAN 
code development. 
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TABLE 1 

SENSOR NUMBER ARRANGEMENT ITERATIONS 
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Subsystem 


B- Positive 
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Figure 21 

Successful Sequential Cube Face Numbering Scheme 
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2 . Algorithm 



The general procedure followed to provide multiple 
layers of redundancy for a force-torque transducer is now 
summarized. It began with a 24 sensor transducer consisting 
of a 2x2 array of FSRs on each of the six inner cube faces. 
The procedures of the previous chapter were utilized to 
determine the normal forces on each sensor resulting from the 
individual application of the three dimensional forces and 
moments. Nominal length dimensions were assumed for d and a 
and the resulting 24x6 matrix was multiplied by a scalar value 
of 24 to produce convenient elements for display purposes. 

The rows of this large matrix were sectioned into four 
invertible 6x6 matrices, arranged sequentially, representing 
individual subsystems capable of being used to resolve the 
force and moment column vector, F, according to Equation 3. 
Each row represented a single sensor, and these subsystems 
were arranged sequentially for FORTRAN programming 
convenience. The ranks of these resulting 6x6 matrices were 
determined to ensure invertibility . 

A sensor was then randomly selected for failure. The 
algorithm had to then determine which subsystem this FSR was 
in and switch to another intact subsystem. This step was 
repeated until no intact subsystems remained. A determination 
was then made as to whether or not there remained enough 
active sensors to conduct a random search for a combination of 
6 remaining sensors to provide a 6x6 invertible matrix for use 
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in resolving F as discussed. If there were at least 6 
remaining sensors, comprised of at least one sensor remaining 
in each the x, y, and z directions, the random search was 
commenced. This search consisted of randomly picking 6 of the 
remaining sensors, again with at least one in each direction. 
The matrix rows associated with these FSRs were then 
sequentially arranged into a 6x6 matrix. The rank of this 
matrix was calculated in order to determine invertiblilty . If 
the rank was six, then the algorithm requested to fail another 
FSR and conduct another search. If the rank was Less than 
six, then a different combination of six of the remaining 
sensors was chosen and the rank was checked. This was 
repeated until a successful 6x6 matrix was produced or a 
sufficient number of iterations had occurred in which case the 
algorithm could be terminated. 

C . FORTRAN PROGRAM 

The flow chart for this algorithm is presented in Figure 
22. The resulting FORTRAN code was written and is presented 
in Appendix B. It was based on the numbering arrangement 
displayed in Figure 21. The values assumed for the length 
dimensions shown Figure 15, d and a, were 4.5 and 1.5, 
respectively. The resulting 24x6 A matrix was multiplied by 
the scalar value 24 in order to produce a simple numerical 
format . 



45 




Figure 22 

Redundancy Analysis Flowchart 
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D. RESULTS 



The resulting FORTRAN program ran successfully. In most 
cases, many combinations were generated in a matter of 
seconds. The slowest part to the process was operator 
interaction and screen display. A sample output is provided 
after the code listing in Appendix B. 

The output shows a listing of the original 24x6 A matrix 
with the option to correct or change any elements. Next the 
sequentially ordered subsystems and their corresponding ranks 
are shown. The program then randomly picks and eliminates a 
sensor. The display shows the sensor chosen and its 
corresponding subsystem. The output then shows which 
subsystems have not yet been disturbed. This is repeated at 
the operators discretion until no more undisturbed subsystems 
remain. At that point the display shows the prograir shifting 
to the mode of randomly picking six of the remaining FSRs and 
generating a new 6x6 subsystem matrix with a rank of six. The 
number of iterations required to generate this subsystem, 
which must meet specific requirements, is atso shown. The 
operator is given the opportunity to allow the program to 
eliminate another FSR and repeat this random generation of a 
useable 6x6 matrix from the remaining FSRs. This process 
continues provided there are at least six sensors remaining 
with at least one in each direction. 
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V. REDUCED ORDER TRANSDUCER THEORY 

A. MOTIVATION 

Many robotic manipulator tasks require only single axis 
moments. An example would be the act of inserting and 
tightening a threaded screw or bolt. The transducer 
configurations considered thus far have been for three forces 
and three moments about each of the x, y, and z axes. In 
order to simplify the transducer design and construction, an 
analogous configuration was developed for a reduced order 
transducer requiring torques about only the z axis and the 
three dimensional forces. 

B. SENSOR CONFIGURATION 
1 . Arrangement 

The concept of a reduced order transducer was 
investigated in the preliminary work which considered each FSR 
capable of detecting only compressive loads [Ref.l:pp. 31-31]. 
However, when the concept of the pre-loaded sensor is 
utilized, the analysis is similar to that used for the full 
order transducer configuration discussed in Chapter III. 

Since it is desired to detect and resolve three forces 
and only one moment, the force and moment column vector, F, 
has only four elements. The reduced order mathematical 
generalization of the problem is therefore: 
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As before, S represents the sensor voltage output and F is the 
force and moment column vector. Also, since it is desired to 
minimize the number of sensors required to resolve the 
components of the 4x1 F vector, n will be four. The goal was 
to choose a configuration which will yield an invertible 
square A matrix. 

Figure 18 showed the modified preliminary transducer 
configuration which was successful in providing an A matrix 
with a full rank. Hence, in order to obtain a reduced order 
transducer, the iterative process of removing sensors and 
their corresponding A matrix rows began with this layout. 

2 . Equations 

As described in Chapter III, the principle of 
superposition was employed to develop the equations describing 
the normal forces on each sensor as a result of the 
application of the three dimensional forces. Also, equations 
were developed to reflect the application of moment about only 
the z axis. The resulting complete equation is similar to 
Equation 22. The only difference is that the two A matrix 
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columns which were a result of the application of M x and My, 
are now gone. The new equation is as follows: 
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3 . Sensor Placement 

First, the values for d and a were assumed to be 5 and 
1, respectively. Then the complete A matrix was put into 
MATLAB™ and the rank was calculated to be four. As expected, 
this indicated a sufficient configuration to resolve the 
components of F. In an effort to reduce the number of sensors 
in the layout, all the sensors on the negative faces of the 
cube were eliminated by setting their associated rows in the 
A matrix to zero. This resulted in a configuration with six 
sensors on the positive face of the cube which yielded an A 
matrix with a rank of four. 
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The layout was then reduced to five sensors by next 
removing one of the two remaining sensors on the positive z 
face. The choice was sensor #11. The associated A matrix 
yielded a full rank of four. Next, sensor #12, the last 
remaining z face sensor, was removed. The resulting A matrix 
was rank deficient indicating a combination incapable of 
resolving F. Hence, this sensor was returned to the 
configuration. 

Sensor #5, ay face FSR, was the next choice. The 
resulting configuration generated an A matrix with a rank of 
four capable of being used to resolve the components of the 
force and moment column vector. For continuity, this sensor 
was returned and sensor #1, a positive x face sensor, was 
removed. This also yielded an A matrix with a rank of four 
indicating another successful combination. Finally, a 
configuration which was nearly a mirror image to the layout 
missing sensor #5 was used generate an A matrix of full rank. 

4 . Results 

After the eight iterations discussed above, two 
useable sensor configuration were determined. These layouts 
consist of the minimum number of pre-loaded sensors, four, 
necessary to resolve the elements of the 4x1 force and moment 
column vector using Equation 13. Figures 23 and 24 display 
the resulting transducer configurations. 
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Figure 23 

Positive Face Pre-Loaded Sensor Reduced Order Transducer 




Figure 24 

Negative Face Pre-Loaded Sensor Reduced Order Transducer 
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VI . PROTOTYPE TRANSDUCER 

A. 4 INCH REDUCED ORDER VERSION 

1. Objective 

It was decided to build a suitable prototype force- 
torque transducer for the reduced order sensor configuration. 
The layout to be used is pictured in Figure 23 and was shown 
theoretically capable of being used to resolve forces in three 
directions and moments about a single axis. After 
calibration, a physical model could be used to qualitatively 
demonstrate the capability of resolving the associated force 
and moment column vector given some arbitrary applied loads to 
the joystick handle. 

2. Design and Construction 
a . Mechanical Design 

Appendix C consists of the detailed machine 
drawings of a 4" force-torque transducer. Figure 25 shows a 
general view of the prototype design. The outer cube was 
constructed using .25'' aluminum plate surrounding a 

plexiglass, smaller solid cube with an attached plexiglass 
joystick handle. The handle was fastened using a pin and 
epoxy. Phillips head 4x1/2" tapping screws were used to 
secure the walls of the outer cube. 
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Figure 25 

Generic Prototype Transducer 



The 4" outside dimension was chosen in order to 
allow for a 3" inner cube using .25" aluminum plate and a .25" 
space between cubes. The need for the 3" inner cube was 
established by the size of the readily available .25" square 
FSR, Interlink Electronics™ product #30-301. Each sensor had 
an overall length of 2". In order to mount these on the inner 
cube and allow enough room for sensor protective padding and 
sensor wiring, the 3" dimension was chosen. 

The four FSRs were configured as shown in Figure 
23, allowing for the detection and resolution of the four 
elements of the force and moment vector according to Equation 
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3. The sensors were mounted to the faces of the inner 
plexiglass cube using Scotch™ 467MP High Performance 
adhesive. A 10/32 threaded set bolt was inserted above each 
sensor site to allow for the application of ar. initial 
compressive force as a pre-load. 

To protect the sensing surface of each FSR, a 
. 5"xl.75" steel shim was fastened to the inside face of the 
outer cube. The sensor side of this shim was initially lined 
with .25" neoprene padding material. This padding was later 
changed to a .30" polyurethane bumpon R to better protect and 
transmit normal forces to the FSR. This alteration provided 
enough padding in the gap between cubes to exert a sufficient 
pre-compressive load to eliminate the set bolts except as a 
method for manually controlling or refining the set point. 
b. Electronic Interface 

The electronic interface shown in Figure 12 was 
used to convert the change in resistance provided by each 
sensor to a suitable voltage change. The sensors were wired 
such that FSR #1 used orange/yellow colored wires, FSR #2 used 
brown/red colored wires, FSR #6 used blue/green colored wires, 
and FSR #12 used gray/ violet colored wires. These were multi- 
strand, 24 gauge wires separated from 16 wire ribbon cable. 
Various feedback resistance values were used in conjunction 
with an LM324 op-amp. This circuit was essentially four of 
the circuits used in the pre-load sensor test apparatus. 
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3. Calibration and Tasting 
a . Apparatus 

A general sketch of the calibration and test 
apparatus is shown in Figure 26. The device was constructed 
using readily available pine l"xlO" material with .125" press 




board used to reinforce the corners as shown. The uprights 
were 8.5" high and the device was 24" wide. The force-torque 
transducer was mounted as shown enabling individual loading in 
each direction of interest. A system consisting of two 
aluminum shafts mounted through ball bearing supports was 
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added to reduce friction between the uprights and load 
supporting strings which pass over the ends. 
b . Technique 

Incremental masses hung from strings laid over the 
upright ends and attached to the joystick handle served to 
load the system in the x and y directions. The joystick cube 
had be rotated 90° for measuring the positive and negative y 
forces after the positive and negative x direction forces were 
applied and measured. 

To measure forces in the positive z direction 
incremental masses were placed in a tray which was mounted to 
the top of the joystick handle. To measure forces in the 
negative z direction, the apparatus was designed to be turned 
over across a gap between two adjacent parallel tables. The 
loads were applied by incrementally placing masses in a tray 
hanging from the joystick handle. 

A bar was bolted to the top of the joystick handle 
parallel to the upright ends in order to measure moments about 
the z axis. Masses were hung in opposite directions from each 
end of this bar using strings passed over the rollers on the 
uprights providing a couple with the bar length as the moment 
arm. The right hand rule was used to determine the directions 
in which to load the bar for the respective posi.tive and 
negative z axis moments. 
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c. Data Acquisition 

Masses ranging from -2.762 kg to +2.762 kg in 
increments of .25 kg were used to calibrate the transducer. 
The moment arm provided by the bar bolted to the top of the 
joystick handle was .0984 meters. The masses were carefully 
used to apply isolated forces in x, y, and z directions as 
well as isolated moments about the z axis. For calibration, 
the OMEGA BENCH™ data acquisition system and APPLE MACINTOSH™ 
microcomputer were configured as a voltmeter to detect the 
voltage response for each sensor as each load in each 
direction of interest was applied. The changes in voltage for 
each change in load were plotted as 16 separate curves. The 
slopes of these curves were used to relate the force and 
moment column vector to the sensor output according to 
Equation 24 as the 16 elements of the 4x4 A matrix. 

With the elements of the A matrix determined, 
Equation 13 could be used to resolve arbitrary applied reduced 
order forces and moments. The OMEGA BENCH™ data acquisition 
icon system was configured as shown in Figure 27 in order to 
perform the matrix multiplication in real time. The first 
column of icons represents the incoming analog signals 
provided by the sensor outputs in response to an applied load 
on the joystick handle. The second column of icons shows 
calculation boxes in which the sensor voltage signals were 
summed with the initial pre-load voltages. This resulted in 
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a change in voltage response to the applied load which could 
be multiplied by the elements of the inverted A matrix. The 
multiplication was done in the third and fourth columns of 
icons by first multiplying each change in sensor voltage 
signal by the inverted A matrix element in the corresponding 
matrix column. Then adding each matrix row of products 
provided each of the four outputs. Each force or moment 
output was displayed in real time using a meter and a chart 
icon. 

d. Results 

Two calibration tests were conducted. The first 
used a feedback resistor, Rq, with a value of 20 kohms and 
the second used 10 kohms for feedback resistance. In both 
tests, it was possible to obtain very linear force and moment 
versus voltage curves, although the curves appeared to be 
somewhat discontinuous at the zero points. This was due to 
the physical variations in first measuring the response to 
applied loads in the positive directions and then unloading 
the axis under consideration and applying the loads in the 
negative direction in order to measure the response. This 
discontinuity was accounted for by determining the slopes 
using four different techniques. The first was slopes 
determined using a least-squares regression line over the 
entire range. The second and third were using the slopes from 
just the positive and negative portions of the curves, 
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Icon Arrangement for Real Time Matrix Multiplication 
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respectively. The fourth was taking the average of the two 
slopes determined for each curve using the second and third 
technique of isolating the respective positive and negative 
ranges . 

Once calibrated using the 20 kohrn feedback 
resistance, the procedure discussed in the previous section, 
and the slopes resulting from the least-squares regression 
line, the system was tested by applying a known force of 10 N 
in the positive x direction. The first time this load was 
applied the response was good. That is the display of forces 
in the y and z directions remained less than about 2 N and the 
moment display registered very little change while the x 
direction force display was about 8.5 N indicating successful 
decoupling of the forces and moments. However, the system did 
not prove to be repeatable. When a force of 10 N was applied 
in the negative x direction the force and moment displays in 
the y and z directions began to show a substantial response 
indicating poor force and moment resolution. The inverted A 
matrix elements resulting from the various slope determination 
techniques were all tried in the data acquisition system in an 
attempt to produce a more repeatable arrangement. None seemed 
to provide a satisfactory result. The matrix elements from 
the technique of determining the slope from only the positive 
load application outputs seemed to qualitatively provide the 
best force and moment resolution. 
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The calibration procedure was repeated using the 10 
kohm feedback resistance with no pre-load set bolts in an 
attempt to refine the data taking technique and overcome this 
non-repeatability. The results were virtually the same. 

B. 2 INCH REDUCED ORDER VER8ION 

1. Objective 

In parallel to building a 4” reduced order force- 
torque transducer, it was desired to design a smaller reduced 
order force-torque transducer to meet the specific payload 
requirements of the PUMA 560 robotic manipulator. This 
required a device with a 2” outside dimension weighing less 
than 5 lbs. In order to do this, a shorter force sensing 
resistor was needed. 

2 . Shortened FSR 
a . Methods 

The smallest FSRs available measured about 1.5" in 
length overall. It was desired to use the same type of FSR 
used in the 4" transducer design in the smaller transducer 
since these provided a good response for a given force 
application. Since the outside dimensions of the transducer 
had to be approximately 2", the inside cube could only be 
about 1.25" to allow a .1875" gap using .1875" aluminum plate. 
Hence, the FSR resistor used could not be any longer than 1" 
to allow enough clearance for wire attached to the terminals. 
Also, the FSRs could be no wider than .5". 
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The .25" square FSRs used in the 4" transducer had 
an overall length of 2". In order to reduce this length, two 
methods were tried which did not use soldering due to the low 
melting temperatures of the materials used to manufacture the 
FSRs. Both methods attempted to leave the sensing pad of the 
device undisturbed and only shorten the tail portion which was 
used to connect to the circuit wiring. The first method 
consisted of separating the outside polymer layers on the tail 
portion of the device. With the silver electrode terminals 
exposed, a 1" section was removed. The remaining .25" portion 
crimped with the wire terminals was inverted and attached to 
the exposed electrodes near the sensing portion with a 
electrically conductive epoxy. After a several attempts an 
adequate bond between the shortened pieces was not achieved. 

The second technique again required splitting the 
polymer covering layers on the tail portion of the FSR. The 
device was trimmed such that only .125" of exposed electrode 
remained as the tail piece. A small, .0625"x.l25" piece of 
the polymer was trimmed away from each side of each exposed 
silver electrode. Small metal terminal pin connectors with 24 
gauge multi-strand wires pre-soldered in position were then 
spread with conductive epoxy and crimped onto the exposed 
electrodes. This provided a strong bond and wire attachment 
for the shortened FSRs. 
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b. Test and Comparison 

To determine the effects of this shortening 
process, a shortened FSR was placed in the pre-load sensor 
test apparatus shown in Figure 11 using the electronic 
interface shown in Figure 12. The test conducted was similar 
to the previous test conducted to prove the pre-load concept. 
A plot of the results is provided in Figure 28. The results 
showed a very linear relationship between the applied force 
and the voltage response when compared to the earlier pre-load 
concept test results shown in Figure 14 . 




Shortened FSR Performance Test Results 
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It was observed that for the shortened test the 
time between incremental loadings was relatively short and 
consistent. During the pre-load concept test, no specific 
care was taken to provide a consistent time interval between 
load applications. This was thought to be important since at 
this point it was noticed that the output slowly crept higher 
as the device was loaded and left undisturbed. 

The test was repeated using a feedback resistance 
of 10 kohms for three different initial pre-load set points. 
A plot of the results is provided in Figure 29. Run #2 and 




Shortened FSR Performance Tests at Various Pre-Loads 

run #3 appear to provide the best set point since they offer 
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the largest change in voltage for a given change in load and 
remain roughly linear. At this point, the shorter FSR appears 
to provide a more linear output than the unaltered sensor. 
This was initially thought to be the result of consciously 
applying the incremental loads and recording the voltage 
response data in a quick consistent fashion in order to 
compensate for the drift phenomenon noticed earlier. 

To test the effects of this data taking technique, 
an unaltered FSR was placed in the apparatus and loaded to 
attain three set points similar to the previous test on the 
shortened FSR. Figure 30 shows the results of this test. The 
response curves are nearly linear indicating that the speed 
and consistency of the data taking technique was important due 
to the drift tendency. The important difference was. that the 
pre-load set points used in the shortened FSR test were not 
repeatable in the unaltered FSR test using a feedback 
resistance of 10 kohms. The higher set point was attained 
using a feedback resistance of 30 kohms while the lower two 
set points were made using a 20 kohm feedback resistance. This 
indicated that the shortening process added about 10 to 20 
kohms of resistance in series with the shortened FSR since the 
unaltered sensor required additional resistance to attain the 
same three set points used earlier and shown in Figure 29. 
The resistance added to the altered FSR is thought to be the 
result of the conductive epoxy used in the fastening process 
of the wire terminals. 
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Unaltered F8R Test at Three Different Pre-Loads 



c. Drift Analysis 
(1) Description 

As discussed in the previous section, the drift 
phenomenon was noticed while applying loads and recording the 
voltage response during testing of the shortened FSR. The 
presence of voltage response drift was confirmed by leaving 
the apparatus loaded and undisturbed for a period of about 3 
hours. In this time the voltage output ranged from 5.988 
volts to 6.290 volts. However, when an arbitrary compressive 
force was suddenly applied and released, the voltage output 
dropped back to approximately 5.9 volts. Therefore, this 
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indicates that the effects of drift do affect the output and 
the response curves whose slopes are used to determine the 
elements of the A matrix. 

(2) Test Conducted 

To more fully understand the effects of drift, 
an arrangement consisting of the four sensors of the 4” 
prototype transducer with a feedback resistance of 10 kohms, 
an unaltered FSR in the pre-load test apparatus with a 
feedback resistance of 30 kohms, an unaltered FSR with no load 
with a feedback resistance of 1 Mohms, a shortened FSR loaded 
with vice grips using minimal padding with a feedback 
resistance of 1 Mohms, and an open channel with 5 volts 
directly applied using no FSR with a feedback resistance of 1 
Mohms was developed. Therefore, this arrangement consisted of 
eight channels, numbered channel 1 to channel 8, representing 
the eight respective sensor setups just delineated. The 
voltage output was plotted over a specific period of time at 
various intervals. 

Two tests were conducted using this eight 
sensor arrangement. The first test was over a period of 
almost 65 hours with 17 output recordings at various time 
intervals for all eight channels with little delay between 
channels. During this test a sudden compressive load was 
applied and released twice and a facility power outage 
occurred. The second test covered a 175 minute period with 
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the output data taken regularly every 10 minutes for each of 
the respective eight channels. 

Another test was conducted using a single 
unaltered .25" square FSR in order to further study the 
effects of the padding used between the load application and 
the FSR active surface. This test consisted of placing the 
FSR between two smooth rigid metal surfaces. One was aluminum 
plate and the other was steel plate material. The data was 
taken first over a period of 276 minutes in 13 intervals. 
These were every 3 minutes for the first 27 minutes and then 
at the 246, 261 and 276 minute points. A second set of 
readings was taken over a 70 minute period in 23 interval. 
These intervals were every minute for the first 10 minutes and 
increased to every 5 minutes for the remainder of the time. 

(3) Results 

The results of the four tests discussed in the 
previous section are provided in Appendix D. The first test 
was conducted over the longest interval. The general 
conclusions from the data plots is that the drift is apparent 
in all channels except channels #6 and #8 which were the 
channels with the unloaded FSR and no FSR, respectively. It 
was also apparent that as the load condition was disturbed, by 
applying a sudden change and quickly returning the original 
loading condition, the voltage value would drop back to the 
initial value and resume drifting upward. This load change 
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was experienced by the apparatus at elapsed times of 10.37 
hours and 57.11 hours. Also, at 59.46 elapsed hours the 
facility experienced a power outage which caused the voltage 
value to drop and resume the upward drift. Throughout all 
these disturbances, the outputs of channels #6 and #8 remained 
nearly constant. 

The second test conducted was similar to the 
first except that the overall time period was shorter and the 
interval between voltage readings was more regular at 10 
minutes each except toward the end of the test. The results 
were as expected, demonstrating the drift in all but channels 
#6 and #8. The only disturbance experienced was an 
intentional system power outage after 155 minutes had elapsed. 
It was at this point that the data acquisition interval was 
decreased to 5 minutes. Except for channels #6 and #8 the 
voltage outputs show the drift upward both after initial load 
application and after the power disturbance. In the resulting 
plots the approach of the effected output curves can be seen 
to approach some apparent asymptotic voltage level. 

In order to understand the effects of the 
padding between the applied load and the FSR sensing surface, 
the last two tests were conducted as discussed in the previous 
section. The results of the tests clearly indicate a voltage 
output gradually rising toward some asymptotic value. No 
disturbances were caused or occurred during these two tests. 
The conclusion was therefore drawn that the type of padding 
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between the load and the FSR sensing surface had no effect on 
the upward drift experienced by the voltage output. 

In the eight channel test, Channel #6 was used 
to show that this drift phenomenon was not caused by the 
circuitry used since no drift was experienced when no FSR was 
used. The conclusion drawn from the results of the Channel #6 
test is that an unloaded FSR does not experience a drift in 
response. Hence, it was shown that only FSRs with applied 
loads experience a circuit output drift upward toward some 
final steady state voltage value. This drift can be stopped 
and started again from an initial voltage output value by 
suddenly disturbing the loading or interrupting the system 
power . 

3. Design and Construction 
a . Mechanical Device 

The design of the 2" force-torque transducer was 
very similar to that of the 4" force-torque transducer. The 
detailed mechanical drawings are provided in Appendix E. The 
primary difference other than the size was the round top and 
bottom to allow enclosure of the device in a cylindrical tube. 
Also, the inner cube was made from aluminum with an aluminum 
joystick handle threaded and inserted into a threaded hole in 
the top. The top of the handle was threaded to allow for 
fastening the bar used to measure the moments about the z 
axis. The size of the inner cube allowed installation of the 
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shortened FSRs exactly on the edges in accordance with the 
reduced order configuration shown in Figure 24. This layout 
was chosen so that pre-load set bolts could be applied to the 
sensor in the z direction at the top of the cube, which is the 
negative z face, in order to allow mounting the transducer on 
the end effector of the PUMA 560 robotic manipulator. 
Polyurethane bumpons R which were .44 ,, x.20 ,, were used between 
cubes in order to apply normal forces to the sensors. Since 
these were slightly larger than the .1875" gap, they provided 
enough pre-compressive load after completely assembling the 
device to eliminate the need for the set bolts. The set bolts 
were therefore only necessary as a means of fine tuning the 
initial set points by adjusting the individual pre-loads on 
each sensor. 

b. Electronic Interface 

The current-to-voltage circuit shown in Figure 12 
was again used for this smaller design. Feedback resistance 
values of 10 and 20 kohms were used. The circuit design was 
essentially identical to that used for the- 4" prototype 
transducer discussed earlier in this chapter. The only 
difference was the wiring color scheme used. FSR #3 used the 
gray/violet colored wires, FSR #4 used the white/black colored 
wires, FSR #7 used the brown/red colored wires, and FSR #10 
used the orange/yellow colored wires. The wires were the same 
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type of 24 gauge multi-strand insulated wires cut from 16 wire 
ribbon cable. 

4. Calibration and Testing 

The technique of calibrating the 2" transducer was 
virtually the same as that of the 4" transducer. The 
apparatus shown in Figure 26 and the OMEGA BENCH™ icon 
arrangement shown in Figure 27 were designed to be used with 
both the 2" and 4” devices in a manner similar to that 
described for the 4” version earlier. Four calibration and 
testing attempts were made. 

The first two calibration tests were conducted using 
a feedback resistance value of 10 kohms. The slopes of the 
resulting 16 least-squared regression lines were calculated 
for each of the two test runs. The respective sets of 16 
slopes were used as the elements of two A matrices. These 
matrices were inverted and used in the calculation icon boxes 
to perform the real time matrix multiplication in accordance 
with Equation 13. The resulting output was similar to that of 
the 4 M transducer when an arbitrary load in -the x, y or z 
direction was applied. The force was not consistently 
resolvable. As the joystick was manipulated, force 
resolution, if attained, was not repeatable. As arbitrary 
moments about the z axis were applied, the same effects were 
observed. The third calibration test was conducted using a 
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feedback resistance of 20 kohms. The results were also not 



repeatable with no clear resolution of forces or moments. 

The fourth test was conducted using no pre-load set 
bolts which had been used earlier to adjust the initial pre- 
load voltage to the nominal starting value used during 
calibration. This test showed very good resolution of 
arbitrarily applied forces in the negative direction, but 
still no repeatable results were attained. 
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VII. DISCUSSION OF RESULTS 



A. DRIFT PHENOMENON 

The presence of drift in output voltage over time was 
demonstrated in the previous chapter. This phenomenon was 
discovered while attempting to calibrate the 2” FSR force- 
torque prototype transducer. The drift in voltage output was 
noticed while loading the device. The output crept upward as 
the load was left undisturbed. The output would drop and 
resume the upward drift from a voltage very near the starting 
voltage as a result of disturbing the load. The load 
disturbances during the tests consisted of suddenly applying 
and removing a larger load. Also, an interruption in power 
would cause the drift to repeat. In all cases the drift in 
voltage appeared to approach some asymptotic value. 

The effects of this drift were not accounted for in the 
initial test of the pre-load concept. The output is shown in 
Figure 14 and does not appear as linear as the output shown in 
Figure 28. The output shown in this later figure was produced 
by attempting to compensate for the drift by using equal 
increments of mass additions over roughly equal intervals of 
time. In an effort to remain consistent, the emphasis was on 
preventing any disturbance which would disrupt the drift while 
taking data. The artificiality imposed by this rigid data 
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taking technique was duplicated during the calibration 
procedures conducted on the prototype transducers. However, 
in actual use of the transducers, it would not be possible to 
accurately account for the drift phenomenon. Sudden changes 
in loading are inherent in the operation of force-torque 
transducers. 

Therefore, it was determined to disregard the effects of 
the drift phenomenon while developing physical models of the 
force-torque transducers based on force sensing resistor 
technology. When a prototype is built which is capable of 
repeatedly resolving the desired forces and moments, the error 
resulting from the drift could be determined. Not only would 
drift effects have been difficult to quantify, during the 
design of the prototypes it was considered superfluous since 
the drift appeared to be less than 10% of the pre-load voltage 
set points. The accuracy of the FSRs used in this application 
were unknown and hence the overall uncertainty of the 
transducer output was not known. At this point in the design, 
accuracy was not the primary concern. The goal was to show 
actual decoupling of the elements of the force and moment 
column vector, F, according to Equation 13. 

B . PROTOTYPE PERFORMANCE 

l. Single FSR Prototype 

The single FSR prototype displayed in Figure 11 using 
the electronic interface shown in Figure 12 was used as a test 
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platform for the pre-load concept discussed in Chapter III and 
for the shortened FSR discussed in Chapter VI. As mentioned 
previously, during the testing for the single pre-loaded FSR 
care was not taken to consistently load the device without 
disturbing the drift inherent in the voltage output. However, 
the device was used to successfully show the ability of the 
pre-loaded FSRs to be used for a voltage signal change, for 
not only compressive loads but also tensile loads. This 
result was significant because it allowed the required number 
of sensors in the transducer design to be halved. This was 
because one pre-loaded sensor was now used to measure forces 
in both the positive and negative directions of the axis on 
which it was mounted. 

The second use of the single FSR prototype was that of 
proving the successful FSR shortening process. A total of 10 
.25" FSRs, Interlink Electronics™ product # 30-301, were 
shortened as discussed in the previous chapter. After each 
FSR was altered, it was mounted in the single FSR pre-load 
device and tested to ensure a adequate response. 

2 . Transducer Prototypes 

The results of the single FSR test device were applied 
to the three dimensional theoretical and physical force-torque 
transducer models. The impact of the successful 
demonstrations of the pre-load concept and the FSR shortening 
process were significant. Also, the effects of the interface 
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between the force application and the FSR sensing surface 
proved to be important. 

First, the pre-load concept was extended to the 
theoretical model of the three dimensional transducer allowing 
the associated A matrix, originally introduced in Equation 2, 
to be square. This was a direct result of halving the number 
of sensors required to detect forces in each direction. It 
was shown for the full order case, which sought the forces in 
each of the x, y, and z directions as well as the moments 
about these three axes that the force and moment column 
vector, F, would have six elements. The elements could be 
resolved using a transducer configuration of only six sensors. 
A similar demonstration was made for the reduced order case 
which sought the forces in the three directions but only the 
moment about a single axis. This simplified the problem, as 
shown in Equation 13, by allowing a square A matrix which, if 
constructed so as to provide a full rank, could be inverted 
and use to solve for F directly. 

Second, the ability to shorten the individual FSRs 
provided the means of building a force-torque transducer of 
smaller dimensions to meet specific robotic manipulator 
payload requirements. The smaller design was analogous to the 
first larger physical three dimensional design except that it 
used the altered FSRs whose performance had been proven on the 
single FSR test device. The FSRs altered in the manner 
discussed in Chapter VI gave the ability to proceed with the 
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design of an adequate force-torque transducer to be used with 
the PUMA 560 robotic manipulator without the time and expense 
of ordering custom FSRs manufactured to meet the specific size 
requirements . 

Finally, in the actual construction, calibration and 
testing of the 2" and 4" prototype transducers, output 
repeatability was a problem. Both prototype devices provided 
an initial response which showed definite decoupling of the 
forces and moments of the F vector according to the reduced 
order theory. However, a consistently repeatable response was 
not obtained regardless of many attempts to accurately 
determine the elements of the A matrix during calibration 
procedures. Ultimately, this lack of repeatability was 
attributed to the tendency of the inner cube to become skewed 
in an orientation with respect to the outer cube which was 
somewhat different than the original orientation. The padding 
used to transfer the load to the sensor was a set of 
polyurethane hemispherical devices which theoretically 
provided a single point of contact on the- FSRs sensing 
surface. In actuality, each pad was compressed a small amount 
providing a relatively small flat spot that remained in 
contact with the sensing surface. In either case, the load 
across the FSR was applied in a relatively small local pattern 
which was free to move across the FSR as the joystick handle 
was moved. Moving the joystick caused the inner cube to move 
with respect to the outer cube. The padding was firmly 
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attached to the outer cube while the FSR sensing surface was 
securely attached to the inner cube. The result was a 
situation in which the application point of the load was free 
to change; therefore, the output of the FSR response was 
reflecting a change in the force application footprint. This 
force application footprint was different than experienced in 
calibration of the device requiring a different A matrix to 
resolve the force and moment vector, F. 

This effect could be proven by utilizing a different 
polyurethane pad above each FSR. These pads should not be 
hemispherical, but cylindrical in shape so that the force 
application footprint on each FSR surface remains constant 
over the relatively small range of motion between the inner 
and outer cubes of the transducer. After eliminating the 
force application footprint as a variable, the device should 
be calibrated as discussed earlier and retested. Upon 
successful and repeatable resolution to the force and moment 
vector, the force-torque transducer accuracy could be 
quantified. The device would then be ready to -be placed in a 
force override rate control system as discussed in the 
preliminary work. 

C . REDUNDANCY ALGORITHM 

The FORTRAN code developed in Chapter IV demonstrates the 
ability of a full order force-torque transducer to be built 
with multiple levels of redundancy. The theoretical 
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development of the redundancy algorithm could also be applied 
to a reduced order transducer requiring moments only about a 
single axis as modelled in Chapter VI. Once proven, the 
physical device could be built with multiple levels of 
redundancy by simply adding more FSRs and a computer interface 
which would conduct the search for useable combinations as the 
need arose. The only limit to the levels of redundancy would 
be imposed by the space limitations dictating the possible 
number of sensors in the physical transducer. 

The full order design used to develop the algorithm 
allows, at a minimum, four layers of redundancy simply by 
using all the available corners of the cube. In the sample 
output in Appendix B, it took 11 sensor failures to exhaust 
the four initial subsystems set up by using 2 4 FSRs and 
requiring six sensors for each subsystem. This left 13 
sensors in an arbitrary configuration from which the search 
for subsequent sets of six sensors would be conducted to find 
useable configurations. In the example output, seven such 
useable systems were found in less than one minute after 
eliminating a sensor each time another set was discovered. 
This example provided for 18 sensor failures out of 24 
possible sensors. The maximum levels of redundancy for this 
configuration is 18 since at least six sensors are required to 
resolve the six elements of the force and moment vector, F. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 



• Force Sensing Resistors offer a practical, robust and 
inexpensive alternative to strain gauge technology for use 
in force-torque transducers. 

• Two, three dimensional reduced order force-torque 
transducers were successfully designed, constructed and 
used to qualitatively demonstrate the ability to resolve 
forces and moments. 

• Three dimensional force-torque transducers are capable of 
multiple levels of redundancy using a simple algorithm to 
search for acceptable configurations given an adequate 
remaining layout of sensors after a sensor fails. 

• The effects of voltage output drift over time will effect 
the overall accuracy of resolving forces and moments using 
FSR based force-torque transducers. 

• Effects of the interface between the force application 
surface and force sensing resistor significantly impact 
the repeatability of the force and moment resolution for 
the force-torque transducer due to position instability 
preventing a consistent force application footprint on 
individual sensors. 



B. RECOMMENDATIONS 

• Conduct a study of different force application surfaces 
and quantify the effects of their footprints on FSR based 
force-torque transducer output repeatability. 

• Use adequate force application surface in existing force- 
torque transducer design to obtain a repeatable resolution 
of the element of the force and moments column vector of 
interest. 

• Develop application software to utilize the resolved 
forces and moments from the FSR based force-torque 
transducer in the force override rate control system. 
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TEST #1 OF SINGLE PRE-LOADED FSR 



(volts) 



RESISTANCE: 4.7 ( kohms ) i - 1 25 



DATA: 

MASS (lbm) VOLTAGE (v) 



m 

I 



V. = 



1.97 


1.49 


1.73 


- 1.47 


-.99 


-.72 


-.51 


-.50 


-.48 


-.24 


.24 


.48 


.51 


.72 


1.01 


1.23 


1.48 


1.73 


2.20 


2.46 


3.18 


3.45 


3.69 


3.93 


4.17 



5.0 


5.65 


5.15 


5.71 


6.26 


6.53 


6.62 


6.63 


6.64 


6.74 


7.01 


7.06 


6.99 


7.11 


7.16 


7.20 


7.26 


7.29 


7.38 


7.45 


7.52 


7.61 


7.62 


7.65 


7.69 



FORCE ( lbf ) : 

mg 
F = — — 
1 32.174 



Response Graph: 
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TEST #2-#6 OF SINGLE PRE-LOADED FSR USING VARIOUS FEEDBACK 
RESISTANCES: 



RESISTANCE: various ( kohms ) i - 1 .. 14 



MASS (g) 

m. = 



DATA: 

VOLTAGE ( 

VI. =V2. = V3. 

i i i 



1149.8 


- 1043.1 


- 936.2 


-724.1 


-494.2 


-49.5 


-.01 


.01 


49.5 


494.2 


724.1 


936.2 


1043.1 


1149.8 



v) 

~V4. =V5 

i i 



5.33 


5.6 


5.87 


6.15 


6.39 


5.36 


5.66 


5.97 


6.27 


6.58 


5.41 


5.76 


6.11 


6.46 


6.80 


5.50 


5.95 


6.40 


6.85 


7.31 


5.56 


6.06 


6.56 


7.06 


7.56 


5.61 


6.17 


6.73 


7.30 


7.86 


5.86 


6.26 


6.87 


7.46 


8.05 


5.81 


6.56 


7.31 


8.07 


8.77 


5.84 


6.62 


7.43 


8.21 


8.77 


5.89 


6.74 


7.59 


8.43 


8.77 


5.90 


6.75 


7.61 


8.44 


8.77 


5.90 


6.74 


7.59 


8.44 


8.77 


5.91 


6.75 


7.59 


8.44 


8.77 


5.92 


6.80 


7.69 


8.59 


8.77 



Where the corresponding 
resistances are (kohms): 

R1 = 2.5 
R2 = 5.0 
R3 =7.5 
R4 = 10.0 
R5 = 12.5 



Force is: 



F. = m.g- 10 ' (N) 



Response Graphs: 




D 2.5 kohms 




D 5.0 kohms 
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Response Graphs cont 





D 7.5 kohms 



° 10 kohms 




F. 

° 12.5 kohms 
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APPENDIX B 
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PROGRAM THESIS 

C ******************************************************************* 

c * 

C * THESIS 

C * NAVAL ENGINEERING 

C * FSR FORCE/TORQUE TRANSDUCER 

C * REDUNDANCY ANALYSIS 

C * 

C * — Charles A. Gunzel-III, LT, USN 

C * 

C ******************************************************************* 

PARAMETER ( N-6 , NRA-N , NCA-N , LDA-N , LDGINV-N , NR- 1 ) 

REAL A (24 , 6) , A1 ( 6 , 6) , A2 ( 6 , 6) , A3 (6 , 6) , A4 ( 6 , 6) , 

+ A5(6,6) , FAC (6, 6) , TOL, AMACH, GINV1 ,GINV2 , GINV3 , 

+ GINV4 , GINV5 , RNO( 1000) , RUNF 

INTEGER COL, ROW , RND , I , J,K, RANK1, RANK2 , RANK3 , RANK4 , 

+ RANKS, FSR(24) , GRP1 , GRP2 , GRP3 , GRP4 , Nl, N2 , 

+ XPOS,XNEG, YPOS, YNEG, ZPOS, ZNEG 
CHARACTER ANS*1 

OPEN (20, FILE- 'THESIS. DAT' , STATUS-' NEW' ) 

LOAD THE MATRIX 

— Use transducer configuration shown in Figure 21. 



DATA (A(1,I) ,1-1,6) /7. 5, 0,0,0, -1,1/ 

DATA (A(2,I) ,1-1,6) /-I . 5 , 0 , 0 , 0 , 1 , -1/ 

DATA (A(3,I) ,1-1,6) /0 , 7 . 5, 0 , 1, 0 , 1/ 

DATA (A(4,I) ,1-1,6) /0 , -1 . 5 , 0 , -1 , 0 , -1/ 

DATA (A(5,I) ,1-1,6) /4 . 5 , 4 . 5 , 3 , 1 , -1 , 0/ 

DATA (A(17,I) ,1-1,6) /4 . 5 , 4 . 5, 3 , -1 , 1, 0/ 

DATA (A(7,I) ,1-1,6) /-7 . 5 , 0 , 0 , 0 , 1 , -1/ 

DATA (A(8,I) ,1-1,6) / 1 . 5 , 0, 0 , 0 , -1 , 1/ 

DATA (A(9,I) ,1-1,6) /0 , -7 . 5, 0 , -1 , 0 , -1/ 

DATA (A(10,I) ,1-1,6) /0, 1.5, 0,1, 0,1/ 

DATA (A(11,I) ,1-1,6) /-4. 5, -4. 5, -3, -1,1,0/ 
DATA (A(23,I) ,1-1,6) /4 . 5 , 4 . 5, -3 , 1 , -1, 0/ 

DATA (A(13,I) , 1-1,6) /7.5, 0,0,0, -1,-1/ 

DATA (A ( 14 , I) ,1-1,6) / -1 . 5, 0 , 0 , 0 , 1 , 1/ 

DATA (A(15,I) ,1-1,6) /0, 7 . 5, 0, 1, 0, -1/ 

DATA (A(16,I) ,1-1,6) /0 , -1 . 5 , 0 , -1, 0 , 1/ 

DATA (A(6,I) ,1-1,6) / -4 . 5 , 4 . 5 , 3 , 1 , 1, 0/ 

DATA (A(18,I) ,1-1,6) /4 . 5, -4 . 5, 3 , -1, -1 , 0/ 

DATA (A(19,I) ,1-1,6) /-7. 5, 0,0, 0, 1, 1/ 

DATA (A(20,I) ,1-1,6) /l. 5, 0,0, 0,-1, -l/ 

DATA (A(21, I) ,1-1,6) /0, -7.5,0, -1, 0, 1/ 

DATA (A(22,I) ,1-1,6) /0, 1. 5, 0, 1, 0, -1/ 

DATA (A(12,I) ,1-1,6) /4 .5,-4 .5, -3,-1, -1, 0/ 
DATA (A(24,I) ,1-1,6) /-4 . 5, 4 .5, -3 , 1, 1, 0/ 

C PRINT THE MATRIX 



15 



WRITE (*,*)'The full matrix is as follows:' 
WRITE (20,*) 'The full matrix is as follows:' 
WRITE (*,*) 

WRITE (20,*) 

WRITE (*,400) ( (A (ROW, COL) , COL-1, 6) , ROW-1, 24) 
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WRITE (20,400) ((A(ROW,COL) , COL-1, 6) , ROW-1,24) 



C VERIFY THE MATRIX 



WRITE (*,*) 

WRITE (20,*) 

WRITE ( * , * ) ' Do you wish to change any elements?' 

WRITE (20,*) ' Do you wish to change any elements?' 

READ (*,300) ANS 
WRITE (20,350) ANS 
IF (ANS.EQ.'Y') THEN 
GO TO 25 
ELSE 

GO TO 35 

25 WRITE (*,*) 'Enter the row and column number of the elemen 1 
WRITE (20,*) 'Enter the row and column number of the elemei 
READ (*,*) ROW, COL 
WRITE (20,*) ROW, COL 
WRITE (*,*) 

WRITE (20,*) 

WRITE (*,100) 

+ 'Enter the corrected value of element (' ,R0W, ' , ' ,CC 

WRITE (20,100) 

+ 'Enter the corrected value of element (' ,R0W, ' , ' ,CC 

READ (*,*) A (ROW, COL) 

WRITE (20,*) A (ROW, COL) 

WRITE (*,*)'Do you wish to change another element?' 

WRITE (20,*) 'Do you wish to change another element?' 

READ (*,300) ANS 
WRITE (20,350) ANS 
IF (ANS.EQ.'Y') THEN 
GO TO 25 
END IF 
END IF 

WRITE (*,*)'Do you wish to view the matrix again?' 

WRITE (20,*) 'Do you wish to view the matrix again?' 

READ (*,300) ANS 
WRITE (20,350) ANS 
IF (ANS.EQ.'Y') THEN 
GO TO 15 
END IP 

35 CONTINUE 



C BREAK FULL MATRIX INTO 4 SUB MATRICES 

DO 37 ROW-1,6 
DO 36 COL-1,6 
A 1 ( ROW , COL ) -A ( ROW , COL) 

A2 (ROW, COL) — A( (ROW+6) ,COL) 

A3 (ROW, COL)— A( (ROW+12) , COL) 

A4 (ROW, COL) -A( (ROW+18) ,COL) 

36 CONTINUE 

37 CONTINUE 

C DISPLAY THE 4 SUB-SYSTEMS 

WRITE (*,*)' The A-POSITIVE face system is:' 
WRITE (20,*) 'The A-POSITIVE face system is:' 
WRITE (*,400) ( (A1 (ROW, COL) , COL-1, 6) , ROW-1, 6) 
WRITE (20,400) ( (AI (ROW, COL) , COL-1, 6) , ROW-1, 6) 
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WRITE (*,*) 

WRITE (20,*) 

WRITE (* , *) 'The A-NEGATIVE face system is:' 

WRITE (20, * ) 'The A-NEGATIVE face system is:' 

WRITE (*,400) ( (A2 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (20,400) ( ( A2 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (*,*) 

WRITE (20,*) 

WRITE ( * , *) 'The B-POSITIVE face system is:' 

WRITE (20,*) 'The B-POSITIVE face system is:' 

WRITE (*,400) ( (A3 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (20,400) ( (A3 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (*,*) 

WRITE (20,*) 

WRITE ( * , * ) ' The B-NEGATIVE face system is:' 

WRITE (20,*) 'The B-NEGATIVE face system is:' 

WRITE (*,400) ( (A4 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (20,400) ( (A4 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (*,*) 

WRITE (20,*) 

C CHECK THE RANK OF EACH OF THESE SUBSYSTEMS 
TOL-lOO * AMACH ( 4 ) 

CALL LSGRR (NRA, NCA, A1 , LDA, TOL, RANK1 , GINVA1 , LDGINV) 
WRITE (*,500) 'The rank of the A-POSITIVE face 
+ system is:',RANKl 

WRITE (20,500) 'The rank of the A-POSITIVE face 
+ system is:',RANKl 

CALL LSGRR (NRA, NCA, A2 , LDA, TOL, RANK2 , GINVA2 , LDGINV) 
WRITE (*,500) 'The rank of the A-NEGATIVE face 
+ system is:',RANK2 

WRITE (20, 500) 'The rank of the A-NEGATIVE face 
+ system is : ' , RANK 2 

CALL LSGRR (NRA, NCA, A3 , LDA, TOL, RANK3 ,GINVA3 , LDGINV) 
WRITE (*,500) 'The rank of the B-POSITIVE face 
+ system is:',RANK3 

WRITE (20,500) 'The rank of the B-POSITIVE face 
+ system is:',RANK3 

CALL LSGRR (NRA, NCA, A4 , LDA, TOL, RANK4 , GINVA4 , LDGINV) 
WRITE (*,500) 'The rank of the B-NEGATIVE face 
+ system is:',RANK4 

WRITE (20,500) 'The rank of the B-NEGATIVE face 
+ system is:',RANK4 

IF ((RANK1.EQ.6) .AND. (RANK2.EQ.6) .AND. 

+ (RANK3.EQ.6) .AND. (RANK4.EQ.6) )THEN 
WRITE (*,*) 

WRITE (20,*) 

WRITE (*,*)'Each of these subsystems is invertible.' 
WRITE (20,*) 'Each of these subsystems is invertible.' 

END IF 
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C INITIALIZE 24 FSR /'S TO ZERO 

DO 39 1-1,24 
FSR ( I ) —0 
39 CONTINUE 

C GENERATE A RANDOM FSR/ FROM 1 TO 24 

GRP 1-0 
GRP2-0 
GRP3-0 
GRP4-0 
XPOS-4 
XNEG-4 
YPOS-4 
YNEG-4 
ZPOS-4 
ZNEG-4 
J-0 

45 J-J+l 
FSR(J)-24*RNUNF() 

C ENSURE FSR/ 0 IS NOT A CHOICE 

IF (FSR(J) . EQ. 0) THEN 
J-J-l 
GO TO 45 
END IF 

C ENSURE THIS IS NOT A REPEAT VALUE 

DO 46 1-1, (J-l) 

IF (FSR(J) . EQ.FSR(I) ) THEN 
J-J-l 
GO TO 45 
END IF 

46 CONTINUE* 

C REPORT THE FSR CHOSEN TO FAIL 

WRITE (*,*) 

WRITE (20,*) 

WRITE (*,500) 'The sensor I will simulate failed is FSR /', 
WRITE (20,500) 'The sensor I will simulate failed is FSR /', 

C TERMINATE IF THERE ARE NOT ENOUGH FSR'S REMAINING 

IF ( J . EQ. 19) THEN 

WRITE (*,*) ' TERMINATION— > Sorry, there are not 
+ enough remaining FSRs. ' 

WRITE(20,*) 'TERMINATION— > Sorry, there are not 
-*■ enough remaining FSRs.' 

GO TO 1000 
END IF 

C CHECK WHICH SUBSYSTEM THIS FAILED FSR IS IN 

IF ((FSR(J) . GE . 1 ) . AND . ( FSR ( J ) . LE . 6 ) ) THEN 
GRP1-1 

WRITE (*,600) 'This sensor is in the A-POSITIVE face system 
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WRITE ( 20, 600) 'This sensor is in the A-POSITIVE face system.' 
ELSE IF ( (FSR(J) .GE.7) .AND. (FSR(J) .LE. 12) ) THEN 
GRP2-1 

WRITE ( * , 600) 'This sensor is in the A-NEGATIVE face system.' 
WRITE(20, 600) 'This sensor i6 in the A-NEGATIVE face system.' 
ELSE IF ( (FSR(J) .GE. 13) .AND. (FSR(J) .LE. 18) ) THEN 
GRP3-1 



WRITE ( * , 600) 'This sensor 
WRITE(20, 600) 'This sensor 
ELSE 
GRP 4-1 

WRITE ( * , 600) 'This sensor 
WRITE(20,600) 'This sensor 
END IF 

C SET THIS FSR SIGNAL TO ZERO 

DO 47 COL-1,6 
A(FSR(J) , COL) —0 . 0 
47 CONTINUE 

C CHECK IF THERE ARE ANY COMPLETE 



is 


in 


the 


B-POSITIVE 


face 


system. ' 


is 


in 


the 


B-POSITIVE 


face 


system . 9 


i8 


in 


the 


B-NEGATIVE 


face 


system . 9 


is 


in 


the 


B-NEGATIVE 


face 


system . 9 



INTACT SUBSETS REMAINING 



WRITE (*,*) 

WRITE (20,*) 

IF (GRP1.EQ.0) THEN 
WRITE (*,*) 'The A-POSITIVE 
WRITE (20,*) 'The A-POSITIVE 
END IF 

IF (GRP2.EQ.0) THEN 
WRITE(*, *) 'The A-NEGATIVE 
WRITE(20,*) 'The A-NEGATIVE 
END IF 

IF (GRP3.EQ.0) THEN 
WRITE (*,*) 'The B-POSITIVE 
WRITE(20,*) 'The B-POSITIVE 
END IF 

IF (GRP4.EQ.0) THEN 
WRITE (*,*) 'The B-NEGATIVE 
WRITE (20,*)' The B-NEGATIVE 
END IF 



face 


system 


is 


still 


active. ' 


face 


system 


i 8 


still 


active. ' 


face 


system 


is 


still 


active. ' 


face 


system 


i8 


still 


active. ' 


face 


system 


is 


still 


active . ' 


face 


system 


is 


still 


active. ' 


face 


system 


is 


still 


active. ' 


face 


system 


is 


still 


active. ' 



I F ( ( GRP 1 . EQ . 0 ) .OR. (GRP2.EQ.0) .OR. (GRP3.EQ.0) .OR. (GRP4.EQ.0) ) THEN 
GO TO 55 
ELSE 

WRITE(*,*) 'None of the 4 initial subsystems remain intact.' 
WRITE (20, *) 'None of the 4 initial subsystems remain intact.' 

GO TO 56 
END IP 



C CHECK IF IT IS DESIRED TO TAKE OUT ANOTHER FSR 

55 WRITE (*,*) 

WRITE (20,*) 

WRITE (*,*) 'Press 1 if you wish to eliminate another FSR?' 
WRITE (20,*) 'Press 1 if you wish to eliminate another FSR?' 
READ (*,*) I 
WRITE (20,*) I 
IF (I.EQ.l) THEN 
WRITE (*,*) 

WRITE (20,*) 
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GO TO 45 
ELSE 

GO TO 1000 
END IF 

56 CONTINUE 

WRITE(*,*)'I will randomly pick 6 of the remaining FSRs t 
+ generate another subsystem.' 

WRITE(20, *) ' I will randomly pick 6 of the remaining FSRs 
+ generate another subsystem. ' 

WRITE ( * , * ) 

WRITE (20, *) 

C CHECK IF THERE REMAINS AT LEAST ONE FSR'S IN EACH DIRECTION 
DO 57 I-1,J 

IF ( (FSR(J) .EQ. 1) .OR. ( FSR ( J) . EQ. 2) 

+ .OR. (FSR(J) .EQ.13) .OR. (FSR(J) .EQ.14)) THEN 
XPOS*XPOS-l 
END IF 

IF ( (FSR(J) .EQ.7) .OR. (FSR(J) .EQ.8) 

+ .OR. (FSR(J) .EQ. 19) .OR. (FSR(J) .EQ.20) ) THEN 
XNEG-XNEG-1 
END IF 

IF ( (FSR(J) . EQ . 3 ) .OR. (FSR(J) .EQ.4) 

+ .OR. (FSR(J) .EQ. 15) .OR. (FSR(J) .EQ.16) ) THEN 
YPOS-YPOS-1 
END IF 

IF ( (FSR(J) .EQ.9) .OR. (FSR(J) .EQ.10) 

+ .OR. (FSR(J) .EQ.21) .OR. (FSR(J) .EQ.22) ) THEN 
YNEG-YNEG-1 
END IF 

IF ( (FSR(J) .EQ.5) .OR. (FSR(J) .EQ.6) 

+ .OR. (FSR(J) .EQ.17) .OR. (FSR(J) .EQ.18) ) THEN 
ZPOS-ZPOS-1 
END IF 

IF ( (FSR(J) . EQ .11) .OR. (FSR(J) .EQ.12) 

+ .OR. (FSR(J) .EQ.23) .OR. (FSR(J) .EQ.24) ) THEN 
ZNEG-ZNEG-1 
END IF 

57 CONTINUE 



IF ( (XPOS.EQ.O) .OR. (XNEG.EQ. 0) 

+ .OR. (YPOS.EQ.O) .OR. (YNEG.EQ.O) 

+ .OR. (ZPOS.EQ.O) .OR. (ZNEG.EQ.O) ) THEN 
WRITE (*,*) 

WRITE (20,*) 

WRITE (*,*) 'TERMINATION— > There are not enough FSRs 
+ remaining in each direction.' 

WRITE (20, *) 'TERMINATION===> There are not enough FSRs 
+ remaining in each direction.' 

GO TO 1000 
END IF 

C TERMINATE IF THERE ARE NOT ENOUGH FSR'S REMAINING 

IF (J.EQ.19) THEN 
WRITE (*,*) 

WRITE (20,*) 

WRITE (*,*) ' TERMINATION===> There are not enough FSRs 

+ remaining to continue.' 
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WRITE (20, *) 'TERMINATION— > There are not enough FSRs 
+ remaining to continue.' 

GO TO 1000 
END IF 

C RANDOMLY PICK 6 OF THE REMAINING FSR'S 

65 Nl-0 
N2-1 

66 CONTINUE 
XPOS-O 
XNEG-0 
YPOS-O 
YNEG-0 
ZPOS-O 
ZNEG-0 

DO 69 K-(J+1) ,24 

67 FSR(K)-24*RNUNF() 

C ENSURE FSR/ 0 IS NOT A CHOICE 

IF (FSR(K) . EQ. 0) THEN 
GO TO 67 
END IF 

C ENSURE THIS IS NOT A REPEAT VALUE 
DO 68 I-1,J 

IF ( ( FSR (I)).EQ.(FSR(K))) THEN 
GO TO 67 
END IF 

68 CONTINUE 

C ENSURE THERE IS ONLY ONE FSR PER FACE 

IF ( (FSR(K) .EQ.l) .OR. (FSR(K) .EQ.2) 

+ .OR. (FSR(K) .EQ. 13) .OR. (FSR(K) .EQ. 14) ) THEN 
XPOS-XPOS+1 
END IF 

IF ( (FSR(K) .EQ.7) .OR. ( FSR(K) .EQ.8) 

+ .OR. (FSR(K) .EQ. 19) .OR. (FSR(K) .EQ.20) ) THEN 
XNEG-XNEG+1 
END IF 

IF ( (FSR(K) .EQ.3) .OR. (FSR(K) .EQ.4) 

+ .OR. (FSR(K) .EQ. 15) .OR. (FSR(K) .EQ. 16) ) THEN 
YPOS-YPOS+1 
END IF 

IF ( (FSR(K) . EQ . 9 ) .OR. (FSR(K) .EQ.10) 

+ .OR. (FSR(K) . EQ . 2 1 ) .OR. (FSR(K) . EQ. 22) ) "THEN 
YNEG-YNEG+1 
END IF 

IF ( (FSR(K) . EQ . 5 ) .OR. (FSR(K) . EQ. 6) 

+ .OR. (FSR(K) .EQ. 17) .OR. (FSR(K) .EQ. 18) ) THEN 
ZPOS-ZPOS+1 
END IF 

IF ( (FSR(K) .EQ.ll) .OR. (FSR(K) .EQ.12) 

+ .OR. (FSR(K) . EQ. 23) .OR. (FSR(K) . EQ. 24 ) ) THEN 
ZNEG=ZNEG+1 
END IF 

IF ( (XPOS1.GT. 1) .OR. (XNEG1.GT. 1) 
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+ .OR. (YPOS1 .GT. 1) .OR. (YNEG1 .GT. 1) 

+ .OR. (ZPOS1.GT. 1) .OR. (ZNEG1.GT. 1) ) THEN 
GO TO 66 

ELSE IF (N1.EQ.6) THEN 
GO TO 70 
END IF 

69 CONTINUE 

C SET UP NEW 6x6 MATRIX 

70 CONTINUE 

DO 72 ROW-(J+l) , K 
DO 71 COL-1,6 

A5( (ROW-J) , COL) — A(FSR(ROW) ,COL) 

71 CONTINUE 

72 CONTINUE 

C CHECK RANK fc PRINT NEW MATRIX IF RANK-6 

CALL LSGRR ( NRA , NCA , A5 , LDA , TOL , RANK5 , GINVA5 , LDGINV) 

IF (RANK5.LT. 6) THEN 
Nl-Nl+1 

IF (N1.EQ.N2* (10**4) )THEN 
N2-N2+1 

WRITE(*,800) 'There have been',Nl,' iterations attempti 
+ find an invertible matrix with the remaining FSRs.' 

WRITE(20, 800) 'There have been',Nl,' iterations attempt, 
+ find an invertible matrix with the remaining FSRs.' 

END IF 
GO TO 66 
END IF 

WRITE (*,*) 

WRITE (20,*) 

WRITE (*,800) 'After' ,N1, ' iterations, the randomly 
+ generated subsystem is:' 

WRITE (20,800) 'After' ,N1, ' iterations, the randomly 
+ generated subsystem is:' 

WRITE (*,400) ( (A5(ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE (20,400) ( (A5(ROW,COL) , COL-1, 6) , ROW-1, 6) 

WRITE (*,*) 

WRITE (20,*) 

WRITE (*,500) 'The rank of this new subsystem is:', RANKS 
WRITE (20,500) 'The rank of this new subsystem is:',RANK5 
WRITE (*,*) 'This subsystem provides an invertible matrix., 
WRITE (20,*) 'This subsystem provides an invertible matrix. 

C CHECK IF IT IS DESIRED TO TAKE OUT ANOTHER FSR 

WRITE (*,*) 

WRITE (20,*) 

WRITE (*,*) 'Press 1 if you wish to eliminate another FSR? 
WRITE (20,*) 'Press 1 if you wish to eliminate another FSR? 
READ (*,*) I 
WRITE (20,*) I 
IF (I.EQ.l) THEN 
J-J+l 
XPOS-4 
XNEG-4 
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YPOS-4 

YNEG-4 

ZPOS-4 

ZNEG-4 

73 FSR(J)-24*RNUNF() 

C ENSURE FSR/ 0 IS NOT A CHOICE 

IF (FSR(J) .EQ.O) THEN 
GO TO 73 
END IF 

C ENSURE THIS IS NOT A REPEAT VALUE 

DO 74 1-1, (J-l) 

IF ( (FSR(I) ) . EQ. (FSR(J) ) )THEN 
GO TO 73 
END IF 

74 CONTINUE 



GO TO 56 
ELSE 

GO TO 1000 
END IF 



C FORMAT LIBRARY 



100 


FORMAT 


(IX, A, 12 ,A, 12, A) 


200 


FORMAT 


(1X,24(1X,F6.3) ) 


300 


FORMAT 


(A) 


350 


FORMAT 


( 10X, A) 


400 


FORMAT 


( (T18,6(F5. 1,3X) ) ) 


500 


FORMAT 


(IX, A, 13) 


600 


FORMAT 


(IX, A, 12, A) 


700 


FORMAT 


(IX, 12) 


800 


FORMAT 


(IX, A, 17, A) 


1000 


CONTINUE 

STOP 

END 
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Sample Output... 



98 



The full matrix is as follows : 



7.5 


0.0 


0.0 


0.0 


- 1.0 


1.0 


- 1.5 


0.0 


0.0 


0.0 


1.0 


- 1.0 


0.0 


7.5 


0.0 


1.0 


0.0 


1.0 


0.0 


- 1.5 


0.0 


- 1.0 


0.0 


- 1.0 


4.5 


4.5 


3.0 


1.0 


- 1.0 


0.0 


- 4.5 


4.5 


3.0 


1.0 


1.0 


0.0 


- 7.5 


0.0 


0.0 


0.0 


1.0 


- 1.0 


1.5 


0.0 


0.0 


0.0 


- 1.0 


1.0 


0.0 


- 7.5 


0.0 


- 1.0 


0.0 


- 1.0 


0.0 


1.5 


0.0 


1.0 


0.0 


1.0 


- 4.5 


- 4.5 


- 3.0 


- 1.0 


1.0 


0.0 


4.5 


- 4.5 


- 3.0 


- 1.0 


- 1.0 


0.0 


7.5 


0.0 


0.0 


0.0 


- 1.0 


- 1.0 


- 1.5 


0.0 


0.0 


0.0 


1.0 


1.0 


0.0 


7.5 


0.0 


1.0 


0.0 


- 1.0 


0.0 


- 1.5 


0.0 


- 1.0 


0.0 


1.0 


4.5 


4.5 


3.0 


- 1.0 


1.0 


0.0 


4.5 


- 4.5 


3.0 


- 1.0 


- 1.0 


0.0 


- 7.5 


0.0 


0.0 


0.0 


1.0 


1.0 


1.5 


0.0 


0.0 


0.0 


- 1.0 


- 1.0 


0.0 


- 7.5 


0.0 


' - 1.0 


0.0 


1.0 


0.0 


1.5 


0.0 


1.0 


0.0 


- 1.0 


4.5 


4.5 


■O • 0 


1.0 


- 1.0 


0.0 


- 4.5 


4.5 


- 3.0 


1.0 


1.0 


0.0 


Do you wish to change 


any elements ? 








N 












The A-POSITIVE face system is : 










7.5 


0.0 


0.0 


0.0 


- 1.0 


1.0 


- 1.5 


0.0 


0.0 


0.0 


1.0 


- 1.0 


0.0 


7.5 


0.0 


1.0 


0.0 


1.0 


0.0 


- 1.5 


0.0 


- 1.0 


0.0 


- 1.0 


4.5 


4.5 


3.0 


1.0 


- 1.0 


0.0 


- 4.5 


4.5 


3.0 


1.0 


1.0 


0.0 


The A-NEGATIVE face system is : 










- 7.5 


0.0 


0.0 


0.0 


1.0 


- 1.0 


1.5 


0.0 


0.0 


0.0 


- 1.0 


1.0 


0.0 


- 7.5 


0.0 


- 1.0 


0.0 


- 1.0 


0.0 


1.5 


0.0 


1.0 


0.0 


1.0 


- 4.5 


- 4.5 


- 3.0 


- 1.0 


1.0 


0.0 


4.5 


- 4.5 


- 3.0 


- 1.0 


- 1.0 


0.0 


The B-POSITIVE face system is : 










7.5 


0.0 


0.0 


0.0 


- 1.0 


- 1.0 


- 1.5 


0.0 


0.0 


0.0 


- 1.0 


1.0 


0.0 


7.5 


0.0 


1.0 


0.0 


- 1.0 


0.0 


- 1.5 


0.0 


- 1.0 


0.0 


1.0 


4.5 


4.5 


3.0 


- 1.0 


1.0 


0.0 


4.5 


- 4.5 


3.0 


- 1.0 


- 1.0 


0.0 


The B-NEGATIVE face system is : 










- 7.5 


0.0 


0.0 


0.0 


1.0 


1.0 


1.5 


0.0 


0.0 


0.0 


- 1.0 


- 1.0 


0.0 


- 7.5 


0.0 


- 1.0 


0.0 


1.0 


0.0 


1.5 


0.0 


1.0 


0.0 


- 1.0 


4.5 


4.5 


- 3.0 


1.0 


- 1.0 


0.0 


- 4.5 


4.5 


- 3.0 


1.0 


1.0 


0.0 
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The rank of the A-POSITIVE face system is: 6 

The rank of the A-NEGATIVE face system is: 6 

The rank of the B-POSITIVE face system is: 6 

The rank of the B-NEGATIVE face system is: 6 

Each of these subsystems is invertible. 

The sensor I will simulate failed is FSR # 4 

This sensor is in the A-POSITIVE face system. 

The A-NEGATIVE face system is still active. 

The B-POSITIVE face system is still active. 

The B-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 

1 



The sensor I will simulate failed is FSR # 3 

This sensor is in the A-POSITIVE face system. 

The A-NEGATIVE face system is still active. 

The B-POSITIVE face system is still active. 

The B-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 

1 



The sensor I will simulate failed is FSR t 18 
This sensor is in the B-POSITIVE face system. 

The A-NEGATIVE face system is still active. 
The B-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR t 22 
This sensor is in the B-NEGATIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR # 16 
This sensor is in the B-POSITIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR # 21 
This sensor is in the B-NEGATIVE face system. 
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The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed ie FSR / 1 

This sensor is in the A-POSITIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR # 20 
This sensor is in the B-NEGATIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR t 19 
This sensor is in the B-NEGATIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR t 2 
This sensor is in the A-POSITIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR / 23 
This sensor is in the B-NEGATIVE face system. 

The A-NEGATIVE face system is still active. 

Press 1 if you wish to eliminate another FSR? 
1 



The sensor I will simulate failed is FSR t 11 
This sensor is in the A-NEGATIVE face system. 

None of the 4 initial subsystems remain intact. 

I will randomly pick 6 of the remaining FSRs to generate another subsystem. 



After 0 iterations, the randomly generated subsystem is: 
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4.5 7.5 

1.5 0.0 

- 4.5 - 7.5 

7.5 4.5 

0.0 7.5 

1.5 1.5 



1.0 1.0 

- 7.5 0.0 

0.0 0.0 

4.5 3.0 

0.0 0.0 

0.0 0.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 



0.0 - 4.5 

- 1.0 0.0 

0.0 1.0 

1.0 - 1.0 

0.0 - 1.0 

0.0 - 1.0 



Press 1 if you wish to eliminate another FSR? 

1 

I will randomly pick 6 of the remaining FSRs to generate another subu 



After 1 iterations, 

- 7.5 

1.5 
- 7.5 

0.0 

- 7.5 

4.5 



the randomly generated 
0.0 - 7.5 0.0 

4.5 4.5 3.0 

- 7.5 0.0 0.0 

- 4.5 4.5 3.0 

- 4.5 4.5 3.0 

1.0 1.0 0.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 



subsystem is: 
-l.o o.o 

-l.o 1.0 

o.o l.o 

l.o l.o 

l.o l.o 

- 4.5 1.0 



Press 1 if you wish to eliminate another FSR? 

1 

1 will randomly pick 6 of the remaining FSRs to generate another subs; 



After 



1 iterations, 
- 4.5 
4.5 
4.5 
- 1.5 
- 1.5 
- 4.5 



the randomly generated 



4.5 - 4.5 - 3.0 

- 1.5 0.0 0.0 

0.0 - 7.5 0.0 

0.0 1.5 0.0 

1.0 1.0 0.0 

4.5 3.0 1.0 



subsystem is: 
- 1.0 - 1.0 

0.0 1.0 

-1.0 0.0 

1.0 0.0 

- 4.5 1.0 

1.0 0.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 

Press 1 if you wish to eliminate another FSR? 

1 

I will randomly pick 6 of the remaining FSRs to generate another subs} 



After 



0 iterations, 
4.5 
- 4.5 
4.5 
- 1.5 
- 4.5 
- 1.5 



the randomly generated 



0.0 

4.5 

4.5 

1.0 

4.5 

0.0 



1.5 

4.5 
4.5 
1.0 
3.0 
0.0 



0.0 

3.0 

3.0 

0.0 

1.0 

0.0 



subsystem is: 

1.0 0.0 

1.0 - 1.0 

1.0 - 1.0 

- 4.5 1.0 

1.0 0.0 

1.0 1.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 



Press 1 if you wish to eliminate another FSR? 
1 
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I will randomly pick 6 of the remaining FSRs to generate another subsystem. 



After 



1 iterations, 
4.5 
-1.5 
4.5 

4.5 

0.0 

1.5 



the randomly generated 



0.0 
4 . 5 
1.0 
4.5 
-7.5 
0.0 



-7.5 

4.5 

1.0 

3.0 

0.0 

0.0 



0.0 

3.0 

0.0 

- 1.0 

- 1.0 

0.0 



subsystem is: 

- 1.0 0.0 

1.0 - 1.0 

-4.5 1.0 

1.0 0.0 

0.0 - 1.0 

- 1.0 1.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 



Press 1 if you wish to eliminate another FSR? 

1 

I will randomly pick 6 of the remaining FSRs to generate another subsystem. 



After 1 iterations, 

0.0 

-4.5 

-4.5 

1.5 

4.5 
-1.5 



the randomly generated 

1.5 0.0 0.0 

1.0 1.0 0.0 

4.5 3.0 * 1.0 

0.0 0.0 0.0 

4.5 3.0 -1.0 

0.0 0.0 0.0 



subsystem is: 

0.0 - 1.0 

-4.5 1.0 

1.0 0.0 

- 1.0 1.0 

1.0 0.0 

1.0 1.0 



generate another subsystem. 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 

Press 1 if you wish to eliminate another FSR? 

1 

I will randomly pick 6 of the remaining FSRs to 



After 



27 iterations, 
-1.5 
0.0 
-1.5 
0.0 

4.5 

1.5 



the randomly generated 

1.0 1.0 0.0 

-7.5 0.0 -1.0 

0.0 0.0 0.0 

1.5 0.0 1.0 

4.5 3.0 -1.0 

0.0 0.0 0.0 



subsystem is: 
-4.5 1.0 

0.0 - 1.0 

1.0 1.0 

0.0 1.0 

1.0 0.0 

- 1.0 1.0 



The rank of this new subsystem is: 6 

This subsystem provides an invertible matrix. 

Press 1 if you wish to eliminate another FSR? 

1 

I will randomly pick 6 of the remaining FSRs to generate another subsystem. 



TERMINATION==»> There are not enough FSRs remaining to continue. 
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APPENDIX C 
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APPENDIX D 
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RESISTANCE: various (kohms) 



I I 17 



DATE/TIME ELAPSED TIME ( hrs ) 



VOLTAGE (v) 



dtg. = t. = 



vl . = v 2 = v 3 . = v 4 . = v 5 . = v 6 . = v 7 . = v 8 . - 

i i i i i i i i 



0 


10.35 


10.37 


11.34 


12.44 


13.39 


14.76 


15.64 


16.67 


23.10 


24.78 


57.08 


57.11 


58.23 


59.46 


62.34 


64.67 



232247 


240908 


240909 


241007 


241113 


241210 


241332 


241425 


241527 


242153 


242334 


260752 


260754 


260901 


261015 


261308 


261527 



6.149 


6090 


6.161 


6.187 


6.087 


-.004 


8.704 


5.059 


6.149 


6.102 


6.176 


6.220 


6.216 


.004 


8.704 


5.059 


5.749 


6.060 


5.975 


6.089 


6.089 


-.004 


8.701 


5.059 


5.963 


6.092 


6.085 


6.170 


6.168 


-.004 


8.709 


5.060 


6.009 


6.113 


6.115 


6.197 


6.193 


.004 


8.714 


5.062 


6.041 


6.128 


6.140 


6.210 


6.211 


-.004 


8.719 


5.062 


6.060 


6.138 


6.164 


6.238 


6.230 


-.004 


8.722 


5.062 


6.074 


6.142 


6.174 


6.248 


6.232 


-.004 


8.724 


5.062 


6.083 


6.147 


6.187 


6.255 


6.240 


-.004 


8.724 


5.062 


6.062 


6.106 


6.182 


6.240 


6.236 


-.004 


8.709 


5.064 


6.060 


6.104 


6.178 


6.234 


6.234 


- .004 


8.704 


5.064 


6.011 


6.050 


6.174 


6.197 


6.253 


-.004 


8.691 


5.060 


5.585 


5.759 


5.745 


5.906 


6.006 


.004 


8.679 


5.060 


5.811 


5.880 


5.929 


6.035 


6.109 


-.004 


8.681 


5.060 


5.852 


5.902 


5.963 


6.060 


6.138 


-.004 


8.684 


5.060 


5.844 


5.887 


5.967 


6.048 


6.121 


-.004 


8.646 


5.057 


5.923 


5.950 


6.031 


6.115 


6.176 


-.004 


8.686 


5.060 




CHI — joystick FSR#1 , lOkohms 

CH2 — joystick FSR/2 , lOkohms 

CH3 — joystick FSR#12 , lOkohms 

CH4 — joystick FSR/6, lOkohms 

CH5 — uncut FSR in beam prototype, 30kohms 

CH6 — uncut FSR w/ no applied load, lMohms 

CH7 — shortened FSR loaded with vice grips, lMohms 

CH8 — no FSR, just 5v directly applied, lMohms 
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CHANNEL 1 CHANNEL 2 





CHANNEL 3 CHANNEL 4 





CHANNEL 5 CHANNEL 6 





CHANNEL 7 CHANNEL 8 
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DRIFT ANALYSIS TEST /2 



RESISTANCE: various (kohms) i =1.2! 



DATE/TIME ELAPSED TIME (min) 

dtg. = t. - vl. 

°i i i 



VOLTAGE (v) 

v2. = v3. = v4. = v5 = v6. = v7. = v8. - 

■ i i i i i i 



282010 


282020 


282030 


282040 


282050 


282100 


282110 


282120 


282130 


282140 


282150 


282200 


282210 


282220 


282230 


282240 


282245 


282250 


282255 


282300 


282305 



0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

155 

160 

165 

170 

175 



5.778 


5.856 


6.031 


6.027 


6.147 


.004 


8.631 


5.058 


5.865 


5.923 


6.089 


6 091 


6.203 


-.004 


8.671 


5.058 


5.872 


5.930 


6.096 


6.100 


6.207 


- .004 


8.676 


5.058 


5.886 


5.937 


6.102 


6.106 


6.213 


-.004 


8.678 


5.059 


5.887 


5.939 


6.104 


6.108 


6.217 


-.004 


8.681 


5.059 


5.891 


5.941 


6.106 


6.113 


6.217 


-.004 


8.681 


5.059 


5.889 


5.941 


6.110 


6.113 


6.221 


-.004 


8.681 


5.059 


5.893 


5.941 


6.110 


6.113 


6.222 


-.004 


8.681 


5.059 


5.895 


5.942 


6.110 


6.115 


6.222 


-.004 


8.683 


5.059 


5.899 


5.944 


6.110 


6.117 


6.224 


-.004 


8.683 


5.059 


5.903 


5.946 


6.110 


6.120 


6.226 


.004 


8.683 


5.059 


5.903 


5946 


6.111 


6.121 


6.226 


-.004 


8.683 


5.059 


5.904 


5.946 


6.115 


6.123 


6.228 


-.004 


8.683 


5.059 


5.904 


5.950 


6.117 


6.125 


6.228 


-.004 


8.683 


5.059 


5.901 


5.946 


6.117 


6.123 


6.228 


-.004 


8.683 


5.059 


5.904 


5.948 


6.117 


6.127 


6.230 


.004 


8.686 


5.059 


5.843 


5.906 


6.087 


6.083 


6.205 


-.004 


8.663 


5.059 


5.884 


5.935 


6.106 


6.104 


6.224 


-.004 


8.68! 


5.059 


5.896 


5.939 


6.110 


6.110 


6.228 


-.004 


8.681 


5.059 


5.893 


5.939 


6.110 


6.110 


6.228 


-.004 


8.683 


5.059 


5.884 


5.935 


6.110 


6.110 


6.228 


-.004 


8.683 


5.059 
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DRIFT ANALYSIS TEST #3; 

RESISTANCE: 1 (Mohms) i = I 13 



DATE/TIME 

dig, r 



031611 


031614 


031617 


031620 


031623 


031626 


031629 


031632 


031635 


031638 


032017 


032032 


032047 



ELAPSED TIME 

t. - 



0 

3 

6 

9 

12 

15 

18 

21 

24 

27 

246 

261 

276 



VOLTAGE (v) 
v7. = 

I 



6.811 


7.015 


7.109 


7.148 


7.168 


7.176 


7.187 


7.238 


7.263 


7.276 


7.498 


7.511 


7.513 



(min) 



CHANNEL 7 
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DRIFT ANALYSIS TEST #4: 

RESISTANCE: 1 (Mohms) i 0 22 



DATE/TIHE ELAPSED TIME (min) VOLTAGE (v) 



dlga ( = 



ta - 



v7a - 

t 



032115 


032116 


032117 


032118 


032119 


032120 


032121 


032122 


032123 


032124 


032125 


032130 


032135 


032140 


032145 


032150 


032155 


032200 


032205 


032210 


032215 


032220 


032225 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 



8.147 


8.214 


8.257 


8.271 


8.295 


8.310 


8.316 


8.324 


8.327 


8.332 


8.334 


8.350 


8 368 


8.388 


8.391 


8.393 


8.401 


8.408 


8.410 


8415 


8418 


8423 


8.425 
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APPENDIX E 
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EXPLODED VIEW 
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